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ABSTRACT 


The  Pauli  and  the  Oberhettinger  asymptotic  expansions  for  the 
diffracted  field  produced  by  the  scattering  of  a  plane  wave  by  a 
wedge  are  compared  analytically  and  numerically,  and  their  range  of 
application  is  extended. 

The  Pauli-Clemmow  method  of  steepest  descents  is  used  to  evaluate 
Sommer f eld's  complex  integral  expression  for  the  total  field  produced 
by  the  scattering  of  a  plane  electromagnetic  wave  by  a  perfectly 
conducting  wedge.  This  method  is  applied  in  a  manner  somewhat  different 
from  that  employed  by  Pauli  and  yields  an  asymptotic  expansion  which 
is  simpler  in  form  and  of  wider  applicability  than  Pauli's  original 
expression.  This  generalized  form  of  Pauli's  expansion  can,  for 
example,  be  applied  to  the  computation  of  the  fields  diffracted  by 
wedges  which  have  exterior  angles  less  than  180  degrees.  It  is  shown 
that  simply  by  rearranging  the  terms  in  this  generalized  Pauli  ex¬ 
pansion  a  generalized  form  of  Oberhettinger 's  asymptotic  expansion 
can  be  produced.  This  generalized  expansion  is  applicable  to  problems 
Involving  wedges  having  exterior  angles  less  than  180  degrees,  and  it  is 
comparable  with  Oberhettinger 's  original  scries  when  the  wedge  angle  is 
greater  than  this  value. 

Several  examples  of  the  scattering  of  a  plane  wave  by  a  wedge  arc 
studied  numerically.  The  superiority  of  the  generalized  Pauli  asymptotic 
expression  over  previously  derived  asymptotic  expressions  is  demonstrated 
in  these  numerical  examples. 

These  asymptotic  expressions  are  used  to  obtain  scalar  diffraction 
coefficients  which  arc  valid  in  the  transition  regions  at  the  shadow 
and  reflection  boundaries. 
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CHAPTER  I 


INTRODUCTION 

This  chapter  contains  a  statement  of  the  problem 
considered  in  the  report,  the  motivation  for  considering 
the  problem,  a  summary  of  the  contents,  and  a  brief  survey 
of  the  literature  on  vredge  diffraction. 

The  Problem 

Let  a  plane  time-harmonic  electromagnetic  wave  of 

arbitrary  polarisation  be  incident  on  a  perfectly  conducting, 

infinite,  two-dincncicnal  wedge;  the  situation  is  depicted 

•  • 

in  Figure  1.  The  field  is  propagating  in  free  space,  and 
its  propagation  vector  k  is  perpendicular  to  the  edge  of 
the  wedge.  When  this  wave  strikes  the  wedge,  it  is  scattered. 
The  sum  cf  the  incident  field  and  scattered  field  is  called 
the  total  field.  This  total  field  can  be  broken  down  into 
two  linearly  polarized  fields,  one  having  its  electric 
vector  parallel  to  the  edge  of  the  wedge  and  the  other 
having  its  magnetic  vector  parallel  to  this  edge.  These 
two  polarization  components  can  be  studied  individually. 

An  exact  mathematical  expression,  called  Gp  in  this 
work,  has  been  developed  to  describe  these  polarization 
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Figure  1,  A  plane  electromagnetic  wave  incident 
on  a  perfectly  conducting  wedge  of  exterior  wedge  angled  • 


3  . 

components,  And  this  can  bo  written  in  a  variety  of  Torino* 

Possibly  the  moot  familiar  of  thcoo  to  the  Infinite  eigen- 

•  • 

function  scries  of  products  of  Bessel  and  trigonometric 
functions*  If  Gp  is  to  be  evaluated  at  points  which  arc 
*  large  number  of  wavelengths  from  the  edge  of  the  wedge, 
then  it  is  convenient  And,  in  fact,  necessary  to  wrxtc  it 
in  the  fora  of  an  asymptotic  expression  in  the  variable  hr, 
Where  k  is  the  propagation  constant  and  r  is  the  perpendicular 
distance  from  the  edge  to  tho  field  point. 

Thorc  aro  two  such  asymptotic  expressions  which  are 

commonly  used  to  compute  Gp  in  such  cases.  One  of  these 

1  J> 

was  derived  by  Pauli}  the  other,  by  Oborhcttingcr.  Pauli's 

expression  is  written  as  a  Fresnel  integral  plus  an  infinite 

series  of  confluent  hypergcometric  functions*  Obcrhettinger 's 

expression  contains  a  Frecncl  integral  and  an  infinite 

•  •  • 

series  of  terms  of  the  form 

sE 

(kr)  .  , 

•  • 

where  N  is  an  integer*  One  purpose  of  thio  work  is  to 

Investigate  the  relation  between  those  two  series*  Are 
thoy,  for  instance,  actually  the  same  ccricn  written  in 
different  formo?  Another  purpose  is  to  determine  which  of 
these  series  in  superior  for  computational  purposes* 
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Each  of  there  series  suffers  from  certain  inadequacies. 

Some  of  these  arc  listed  below. 

•  .  • 

ft)  Both  series  arc  restricted  to  applications  where  the 
exterior  wedge  angler'  is  greater  than  180  degrees. 
b>  The  Pauli  series  cannot  bo  applied  to  all  situations 
even  when  of  is  greater  than  180  degrees.  For  example, 
if  oi  is  close  to  but  greater  than  l80  degrees  and  the 
propagation  vector  1:  in  nearly  parallel  to  one  of  the 
faces  of  the  wedge,  the  Pauli  aeries  will  not  describe 
Op  accurately  in  the  vicinity  of  the  other  face, 
c)  The  Pauli  series  is  written  in  terms  of  the  relatively 
unfamiliar  confluent  hyperge one trie  functions. 

In  this  work,  these  series  are  rc-cxoTinod ,  and  these 
Inadequacies  aro  in  part  removed. 

The  Motivation 

The*  geometrical  theory  of  diffraction  provides  a 
simple  method  for  calculating  the  fields  produced  by  the 
scattering  of  high  frequency  electromagnetic  waves  by 
perfectly  conducting  bodies  of  quite  general  shape.  ThiG 
theory  is  particularly  useful  for  describing  the  scattering 

by  edges  of  bodies,  such  as  the  curved  edge  of  a  circular 

4  5 

disk,  the  curv*xl  edge  of  a  parabolic  antenna,  the  straight 

edge  of  a  horn  antenna,^  or  the  straight  edge  of  a  polygonal 

9 

cylinder. 
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In  the  geometrical  theory  of  diffraction,  the  deccrip- 
tion  of  diffraction  from  edges  is  treated  in  the  following 
way.  As  in  the  theory  of  geometrical  optics,  the  field 
incident  on  the  edge  is  assumed  to  propagate  along  ray 
paths.  Upon  striking  an  edge,  the  incident  ray  is 
scattered,  giving  rice  to  a.  family  of  diffracted  rays. 

The  direction  of  propagation  of  these  diffracted  rays  is 
determined  by  the  generalized  Fermat's  principal.  The 
distribution  of  the  incident  energy  among  these  diffracted 
rays  is  described  by  a  diffraction  coefficient.  The 
diffraction  of  the  incident  rays  is  a  local  phenomenon, 
so  the  diffraction  coefficient  depends  on  the  geometry 
of  the  edge  in  the  region  surrounding  the  point  of  diffrac¬ 
tion  and  on  the  local  polarization  state  of  the  incident 
field. 

The  exact  form  of  the  diffraction  coefficient  must  be 
derived  from  the  rigorous  solution  of  a  "canonical"  problem. 
The  canonical  problem,  which  yields  the  diffraction 
coefficient  for  an  edge,  is  the  diffraction  of  a  linearly 

polarized  plane  wave  by  a  tv/o-dimensional,  infinite  v/edge. 

8 

The  diffraction  coefficient  is  found  from  Sommer fe Id ' s 
asymptotic,  high  frequency  solution  to  this  problem. 

This  edge  diffraction  coefficient  adequately  describes 
the  diffracted  field  at  points  well  removed  from  the  point 


of  diffraction  and  transition  regions,  that  is,  shadow  and 
reflection  boundaries.  In  these  transition  regions,  however 
the  diffraction  coefficient  fails  to  describe  the  field  and, 
the  geometrical  theory  must  be  supplemented  by  a  more  cophis 
ticated  equation.  The  equations  normally  employed  arc 
the  asymptotic  expansions  of  the  wedge  diffracted  field 
given  by  Pauli  and  by  Oberhettinger . 

When  applying  these  formulas  to  practical  geomc ti'ical 

theory  of  diffraction  computations  a  number  of  questions 

naturally  arise.  Which  formula  is  superior  from  the 

computational  point  of  view,  or  which  asymptotic  expansion 

will  provide  the  best  description  of  the  field  in  the 

transition  region  v/ith  the  minimum  amount  of  computational 

offort?  Y/hat  is  the  smallest  value  of  kr  for  which  these 

formulas  are  useful?  Since  the  expressions  are  asymptotic 

series,  what  is  the  effect  of  higher  order  terms?  This 

study  was  motivated  by  practical  computational  questions 

• 

of  this  nature  as  v/ell  as  by  the  desire  to  extend  the 
useful  range  of  Pauli  and  Oberhettinger  asymptotic 
expressions  and  to  determine  the  relation  between  them. 

Summary 

Tho  starting  point  of  this  work  is  the  eigenfunction 
series  form  of  the  Green’s  function  for  a  line  source 
radiuting  in  the  vicinity  of  an  infinite,  two-dimensional 


? 


wedge*  After  a  brief  discussion  of  thic  expression,  it  is 

assumed  that  the  line  source  recodes  to  infinity,  and  this 

•  • 

solution  5s  reduced  to  that  for  the  scattering  of  a  plane 
wave.  The  resulting  infinite  scries  is  then  summed  in 
closed  form  and  written  as  a  complex  contour  integral. 

This  integral  is  identical  with  the  an sat a  first  deduced 

g 

by  S  omtscrfcld.  This  complex,  integral  expression  for  Gp 

is  then  rewritten  in  a  slightly  different  but  more  convenient 
form,  and  distribution  of  the  saddle  points  and  poles  of 
the  integrand  of  the  resulting  expression  is  studied. 

The  integral  expression  is  then  evaluated  by  the  ordinary 
Bethod  of  steepest  descents,  and  then  by  a  nodiiicaticn  of 
this  method.  This  modified  method  of  steepest  descents  was 
first  used  by  Pauli  in  the  derivation  of  his  asymptotic 
scries  for  Gp,  The  method  ic  used  here  in  a  somewhat  differ¬ 
ent  manner  to  yield  an  expression  for  Gp  which  is  both 
simpler  and  more  general  than  Pauli’s  expression.  The 
expression  derived  in  this  work  consists  of  an  infinite 
scries  of  functions.  Each  function  is  the  sum  of  a  Fresnel 
integral  and  a  finite  number  of  terms  of  the  form 
-  W 

(Hr)*  . 

It  is  next  demonstrated  that  for  certain  cases  the 
expression  for  Gp  developed  5n  this  work  can  be  reduced 
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to  the  one  derived  by  Pauli.  The  cases  for  which  this  reduction  is 
not  possible  are  just  those  to  which  the  Pauli  expression  does  not 
apply. 

The  relation  between  the  expression  for  Gp  developed  in  this 
work  and  that  developed  by  Oberhettinger  is  investigated.  It  is 
shown  that  for  cases  where  oL  is  greater  than  180  degrees  they  are 
the  same  series  written  in  a  different  form.  In  Chapter  IV  several 
examples  of  scattering  of  a  plane  wave  by  a  wedge  are  studied  numeri¬ 
cally,  and  some  conclusions  regarding  the  relative  merit  of  the  vari¬ 
ous  solutions  presented  in  Chapter  III  are  drawn  from  the  results. 

In  Chapter  V  the  asymptotic  solutions  are  used  to  obtain  expressions 
for  the  diffraction  coefficients. 


The  problem  of  diffraction  of  electromagnetic  waves  by  a  wedge 

has  been  the  subject  of  numerous  theoretical  and  experimental  studies. 

An  extensive  list  of  many  of  the  papers  is  given  in  the  bibliography. 

9 

Sommerfeld  was  the  first  to  develop  a  rigorous  solution  for  the  dif¬ 
fraction  of  a  plane  wave  by  a  perfectly  conducting  half-plane.  He  used 
an  extension  of  image  theory  to  deduce  an  integral  solution  for  this 
problem  and  showed  that  his  solution  reduced  to  the  Fresnel  integral.  A 
similar  solution  was  derived  by  Carslaw,^  who  used  image  theory,  and  by 
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Poincare7,*'*  who  reduced  the  diffraction  problem  to  one  of 

IP  ix 

oplving  an  integral  equation.  Macdonald,  *  using  the 

classical  methods  of  separation  of  variables,  was  the  first 

to  obtain  the  complete  solution  to  the  problem  of  diffraction 

of  plane,  cylindrical,  and  spherical  waves  by  a  perfectly 

conducting  wedge  of  arbitrary  angle.  Macdonald's  solution 

is  in  the  form  of  an  infinite  cerics  but  has  been  written 

in  integral  form  by  means  of  an  integral  transformation 

relating  the  ordinary  and  modified  Bessel  functions. 

8 

Sommcrfcld  employed  his  extension  of  the  method  of  images 

to  derive  an  integral  expression  for  the  diffraction  of 

a  plane  wave  by  a  perfectly  conducting  wedge  of  arbitrary 

angle.  Carslaw* f  and  Bromwich*'*  also  constructed  solutions 

16 

to  this  problem  using  the  method  of  images.  Obcrhettinger 
used  an  integral  transform  method  to  derive  another  integral 
expression  for  the  diffraction  of  plane,  cylindrical,  and 
spherical  waves  by  a  perfectly  conducting  wedge  of  arbitrary 
angle.  Kontorovich  and  Lebedev  (  and  Karp  have  developed 
the  solution  for  a  half-plane  ucing  transform  methods. 

The  Wicncr-Hopf  method  has  been  applied  to  this  problem 
by  Harrington,*^  Copcon,^  and  Magnus.^*’  Recently  Nomura,^ 

pi i  • 

Senior,  J  and  HoinG  have  studied  the  diffraction  of  a  di¬ 
polo  field  by  a  perfectly  conducting  half-plane.  The 
solution  for  the  diffraction  of  pulses  by  wedges  has  been 
examined  by  Otorhct linger,^*  FrcdcUandcr,*^  and  others.^1  *** 
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The  problems  of  diffraction  by  imperfectly  conducting 

wedges,  dielectric  wedges,  and  wedges  having  mixed  boundary 

conditions  have  been  treated  in  recent  year c.^*  ^  Senior 

hac  developed  an  exact  solution  for  the  imperfectly  conducting 
42  43 

half-plane  and  right-angled  wedge.  Approximate  solutions 

44 

to  this  problem  have  been  found  by  Jones  and  Fidduck, 

Fclocn*^  Williams,*1^*  ^  and  Malyuzinec.^  Diffraction 

froin  unidirectional  conducting  half-planes  has  been  studied 

cq  51  52  55 

by  Radlow,^  Hurd,  and  Scchndri,  1  y  and  the  problems 

of  diffraction  by  wedges  immersed  in  conducting  and  anisotropic 

5/1 .  57  55 

media  have  been  investigated  by  Dmitriev,  1  Jull, 

Williams,^  and  Seshadri.^  The  diffraction  of  surface 

vra ves  by  right-angled  wedges  was  treated  by  Karp  and 

Knra.1,^®"^  Chu  and  Kouyoumjinn arid  others. 

Diffraction  of  waves  by  v/edgee  with  various  kinds  of  exotic 

71-75 

boundary  impcdanceo  have  been  examined  by  Karp, 

Felccn,^**^  and  others.®®"^ 

Jrtic  results  of  come  of  the  theoretical  investigations 

of  wedge  diffraction  have  been  applied  to  the  calculation 

91  92 

of  horn  and  surface  wave  antenna  patterns,  *  .to  studies 

93 

Of  scattering  from  terrain  features,  and  objects  possessing 
7  94  95 

edges,  •  *  and  to  the  problems  of  diffraction  from  thick 

ecrconc,^  irregular  edges, 97”^  and  various  double  wedge 

99  100 

configurations.  *  A  number  of  experimental  studies 

have  been  made  to  verify  the  predictions  of  wedge  diffraction 
theory. 101006 


n 


Asymptotic  forms  of  the  solution  for  the  diffraction 

of  a  piano  wave  by  a  perfectly  conducting  wedge  have  been 

o  1  107 

derived  and  studied  by  Sommer fold*  Pauli,  Clemmow, 

Yan  dor  lYaerdcn,^0^  and  Oberhettinger Soraaerfeld^ 

attempted  to  derive  an  asymptotic  expression  for  his  integral 

solution  of  thi6  problem.  However,  his  result  failed  to 

predict  the  correct  value  for  the  field  near  the  shadow  and 

reflection  boundaries.  This  failure  was  due  to  the  fact 

that  he  used  the  conventional  method  of  steepest  descents 

109 

developed  by  Debye.  This  method  is  not  applicable  to 

coses  in  which  a  pole  lies  near  the  saddle  point  over  which 
the  integration  is  taken,  a  situation  which  occurs  when  the 
diffracted  field  is  evaluated  near  a  shadow  or  reflection 
boundary.  Pauli^  modified  the  method  of  steepest  descents 
and  derived  a  suitable  asymptotic  solution  for  the  problem. 
Ott^^  examined  Pauli's  method  and  simplified  it.  Cleniwow*^ 
studied  Pauli's  method  in  detail  and  found  it  yields  a  partial 
asymptotic  expansion.  Clemmow  suggested  that  such  an 
expansion  can  always  be  rearranged  in  the  form  of  an  anymptoti 
expansion  in  inverse  powers  of  an  appropriate  variable. 
Oberhettinger  applied  Watson's  lemma  to  his  integral 
expression  for  the  diffraction  of  a  piano  wave  by  a  wedge 

1  o  3 

to  produce  a  new  asymptotic  expansion.  Van  dor  Waerden 
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and  Obcrhet tinker  showed  that  Obcrhettingoi  *s  expansion 
could  be  derived  by  an  application  of  Van  dcr  Y/aerdcn's 
modified  method  of  steepest  descents. 
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CHAPTER  II 


SERIES  AND  INTEGRAL  FORMS  OF  THE  SOLUTION 

In  thin  chapter  the  in finite-series  form  of  the 
Grcen’c  function  for  the  diffraction  problem  of  a  line 
courcc  radiating  in  the  presence  of  a  two-dimensional 
ircdge  of  arbitrary  angle  is  introduced,  and  the  infinite- 
series  form  of  the  solution  for  the  diffraction  of  a  plane 
wave  by  a  wedge  is  derived  from  this.  This  series  is  then 
transformed  to  a  complex  integral  representation. 

Scries  Form  of  the  Solution 

Consider  a  two-dimensional  wedge  and  a  line  source 
to  be  situated  in  space  as  shown  m  Figure  ?.,  The  faces 
of  the  wedge  are  formed  by  two  corni-infinite  planes 
intersecting  on  the  z-axis  of  the  cylindrical  coordinate 
oyctern.  One  plane  is  located  ut  s  0;  the  other,  at  r.c 
The  infinitely-long  line  source  in  parallel  to  the  edge  of 
the  wedge,  and  itG  position  is  described  by  the  coordinates 
(r'lpi1).  The  typical  field  point  is  denoted  by 

The  lino  source  ic  assumed  to  have  Unit  strength  and  time 

.  <10 1 

Dependence  of  the  form  c‘  • 
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Figure  2. 
✓  . 


A  line  eourcc  and  diffracting  wedge. 


•- .  -  •-  .  .  - 


The  Green's  function  for  thio  radiating  system  can  bo 

l6 


written 


co 

» ><kr)  [coa  t  «(,/j  * 


w  c  0  n 


♦  cos 


[?<**♦  *»•>)] 


if  r  >  r*.  For  the  caoe  r'<  r,  r  and  r'  aro  interchanged* 

In  thio  exprcscicn  Cm  in  the  Neumann  number  which  is  equal 

to  one  if  n  is  i»cro;  otherwise  it  is  equal  to  two.  J^(kr) 

n 

and  IP  (kr)  represent  Bessel  and  Hankcl  functions  respectively. 


The  variable  n  describes  the  wedge  angle  and  io  related  to 
by  the  following  equation: 

*  «  nir  .  (2) 

k  io  the  propagation  constant.  The  plus  sign  is  used  between 
the  two  cosine  terms  if  the  boundary  condition  is  of  the 
homogeneous  Neumann  type  ( *  0  on  both  faces  of  the 
wedge).  For  the  homogeneous  Dirichlct  boundary  condition 
*  0  on  both  faces  of  the  wedge ),  the  ninus  sign  in  used. 
This  convergent  eerier,  is  on  exact  solution  to  the 
time -harmonic ,  inhomogeneous  wave  equation  for  the  problem 
of  a  radiating  line  source  and  wedge  embedded  in  a  linear, 
isotropic,  honogcncouo,  lossless  medium.  It.  satisfies  the 
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homogeneous  Neumann  or  Dirichlot  boundary  conditions, 

11/| 

the  radiation  condition  and  the  Mcixncr  edge 

condition?*'*"*^  This  series  converges  for  all  values 

of  **i  r'i  f * »  and  j'1,  and  is  valid  for  all  values  of  n* 

The  expression  describes  the  total  radiation  field 

created  by  the  direct  radiation  from  the  lino  source  and 

that  scattered  from  the  v/edge.  If  the  line  source  is  an 

electric  current  of  strength  I,  then  G  represents  the 

electric  field  vector,  that  is,  E  s»ji£K7.G;  and  the 

boundary  condition  G  e  0  applies.  If  the  line  source  is 

a  magnetic  current  of  strength  M,  then  G  represents  the 

magnetic  field  vector,  that  is,  li^  =  - jt&MG }  and  the 

d  G 

boundary  condition  ~r  =.  0  applies. 

d  p 

In  many  cases  it  is  necessary  to  determine  the  total 
radiation  field  when  the  field  point  is  far  removed  from 
the  vertex  of  the  v.edge.  In  such  cases  Equation  1  can  be 
cinplified  by  replacing  the  Hankcl  functions  by  the  first 
term  t f  their  asymptotic  expansion,  that  is,  by  the  relation 


TT  ra  IT  \ 
-  H  "  n  2  ) 


This  substitution  reduces  G  to  the  form 


G 


-j(kr  -  £) 


oo 


m  & 


6  J  (ltr')c 
m  m 

mm 

0  n 


<  m  IT 
i  "  2 


n 


5*hic  substitution  is  strictly  valid  only  for  those  terms  of 
the  series  for  which  kr  is  largo  compared  to  the  order  of 
the  Jlankel  function.  Therefore,  when  making  this  substitu¬ 
tion  an  additional  assumption  is  required.  This  is  that 
hr'  is  snail  enough  compared  to  kr  so  that  the  terms  of 

• 

the  scries  which  do  not  satisfy  the  requirement  kr  ^  — 
arc  negligibly  small  and,  for  practical  purposes,  contribute 
nothing  to  C. 

If  the  line  source  rather  than  the  field  point  is  far 
removed  from  the  vertex  of  the  v/edge,  G  can  be  written 


rhcrc 


(5) 


i  coc. 


ft 


(6) 


Cp  is  the  series  form  of  the  Green's  function  describing 
the  total  field  created  by  the  scattering  of  a  plane  wave 
by  a  wedge.  It  ic  an  exact  solution  to  this  problem  in 


The  remainder  of 


the-  Dense  described  on  pages  15  and  16. 
this  work  concerns  the  evaluation  of  Gp. 

Equation  6  converges  rapidly  for  small  values  of  kr. 

If  kr  is  leer,  than  1.0,  less  than  15  terms  of  the  series 
ore  required  to  compute  a  value  of  Gp  which  is  accurate 
to  5  significant  figures.  When  kr  is  large,  the  series 
still  converges  but  very  slowly.  If  kr  u  10,  *10  terms 
ore  required  to  achieve  5  significant  figure  accuracy  in 
Cp.  In  view  of  this  slow  convergence  for  large  values  of 
kr  the  expression  for  Gp  must  be  cast  in  another  form  if 
it  is  to  bo  useful  for  computational  purposes.  An  asymp¬ 
totic  expansion  of  Gp  in  inverse  powers  of  kr  is  such  a 
useful  form.  In  order  to  derive  an  asymptotic  expression 
for  Gp  by  the  standard  method  of  steepest  descents  it  must 
first  be  transformed  into  an  integral  or  integrals 
of  the  type 

(  F(Z)ckrf(z)d*  (7) 

c 

In  the  next  section  a  method  for  making  this  transformation 
Jo  described. 

_  •  * 

Intcrrol  Form  cj*.  tli_£  Solution 

Sommcrfcld  deduced,  by  an  extension  of  the  method  of 
images,  a  complex  integral  form  for  Gp  and  derived  the 
eigenfunction  series  given  by  Equation  6  from  this  integral. 
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In  this  section  the  inverse  derivation  is  presented*  This 
procedure  is  more  systematic  since  one  is  not  required  to 
start  with  a  solution  arrived  at  by  intuition.  First  one 
writes  the  Bessel  functions  of  Equation  6  in  integral  form. 
Next  the  order  of  integration  and  summation  is  interchanged, 
end,  finally,  the  resulting  sums  are  written  in  closed  form. 

Gp  will  first  be  written  as  the  sum  of  tv/o  terms  as 
shown  below: 


Gp  «  Kkr,  -fS‘ ,  n)  +  Kkr,  0  +  ' ,  n)  ,  (8) 

where 


CO 


Kkr,  ±  0*,  n)  =  ~ 


.  mlT 
J  «  D 


n  Y  €  r»Jn(kr)C  "  **  1  > 

Ls  n 


(9) 


m 


In  the  analysis  to  follow  both  of  these  terms  will  be 
denoted  by  I(kr,  p  ,  n),  where  |3  =  ^  , 

If  cos  «  p  is  expressed  in  complex  form,  Equation  9 
can  be  written 


(10) 
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The  BcgscI  functions  in  this  expression  can  be  expressed 

in  terns  of  contour  integrals  in  the  complex  z-plnno  as 
•  • 

follows: 


J  (fcr) 
n 

mm 

n 


jkr  cos  z 
c 


(z 


or 


J.Ckr) 

13 

n 


jkr  cos  z  ~  (z  +  vf) 
e  dz 


(11) 


(12) 


The  contours  c  and  c'  arc  pictured  in  Figure  The 
cross-hatched  strips  in  the  figure  represent  the  regions 
in  which  the  integrands  vanich  when  |lia»T*-00  •  Mow,  if 
tho  integral  over  c  iG  substituted  into  the  first  simulation 
in  Equation  10,  the  integral  over  c'  is  substituted  into 
the  second,  and  if  the  order  of  integration  and  sucks 1 5  on 
is  interchanged  in  each  tern,  the  following  results: 

00 


Kkr.p 


jkr  cos  z 
c 


i  B  <(U  z) 


dz  + 


1  ( 

Sff  ) 


Jkr  cos  z  f  -  J  J  ^  *  "> 

Z< 


dz  (13) 


n  k  0 


z 


1 


Z- PLANE 


7,-planc  and  tv/o  possible  integratis: 
function. 
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'/he  two  emanations  can  be  reduced  to  the  closed  forms 


OD 

2 


i  »  <f  +  *> 


-i 


O  »  1 


1  -  C 


-  3 


St*.*. 


n 


<1*0 


and 

CO 

£■ 

Be  0 


•  J  j(^  +  s) 


1  -  e 


-  3 


.GJL2, 


n 


(15) 


Equation  lj>  can,  therefore,  be  written 


I(kr,£>  ,  n)  .  ^  ( 


(«'  -  <=) 


1  ••  e 


-  3 


ft  h  7. 

n 


ccs  "dr.  (If) 


The  negative  sign  before  the  c  indicates  that  this  integration 

f 

contour  is  to  be  traversed  in  the  direction  opposite  to  that 
indicated  in  Figure  J.  Substitution  of  Equation  16  into  8 

yields 


a  Z 
J  n 


0  -  —i— 

P  21<n 


<  5.  ..  <  SLzX'l. 

v  9%  ^  n 


,jl:r  eon  »  d. 


(c'  -  c) 


n 


-  e 


4  JL 

1  2lTn 


3  * 
J  n 


(c*  -  c) 


3  - 

J  n 


-  c 


J  II 


(jkr  cos  8  da 


(17) 
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Thie  is  the  expression  deduced  by  Sommcrfcld. 

Equation  16  can  be  put  in  a  fora  i/ioro  convenient  for 
later  use  by  noting  that 


3 

J  2n 

_ _ e 

<  iLtS  IT1EIII 

J  2n  J  2n 

o  —  c 


(18) 


Substituting  this  into  Equation  16  and  making  use  of  the 
fact  that 


|  (  cjkp  cos  ziz  s  1  j^(kr)  .  |  Jo():r)  =  0  ,  (19) 

(c1  -  c) 
yields 


(20) 


G„  lc  now  cxproeccd  in  the  integral  form  opccified  by 

V 

Equation  ?•  In  the  next  chapter  two  modifications  of  the 
method  of  steepest  descents  will  be  utilized  to  evaluate 
this  integral  asymptotically. 


CHAPTER  III 


ASYMPTOTIC  EVALUATION  OF  TilE  INTEGRAL  FORM 
OF  THE  SOLUTION 

In  the  previous  chapter  the  eigenfunction  series 
solution  for  the  diffraction  of  a  plane  wave  by  a 
wedge  of  arbitrary  angle  was  expressed  as  the  sum  of 
two  integrals,  each  of  the  form 

I(kr,^,  n)  =  ]~jrr  (  F(z)ekrf^dz  (21) 

(c*  -  c) 

where 

f(z)  =  j  cos  2  (22) 

and 

TM  =  (2j» 

✓ 

In  this  chapter  these  integrals  arc  evaluated  for  the  case 
of  )  arge  hr  by  means  of  the  method  of  steepest  descents 
and  the  I’auli-Clemmow  modification  of  this  method.  In  the 
first  section  the  saddle  points  of  f(z)  arc  located,  and 
the  integration  contour  c  -  c 1  is  closed  by  moans  of  steepest 
descent  paths  through  two  of  these  saddler..  Next,  the  poles 
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of  F(s)  are  located,  their  contribution  to  Gp  is  determined, 
and  their  behavior  for  various  values  of  p  and  n  is  described 
In  the  next  section  the  integrals  which  remain  in  I(kr,f,  n) 
are  transformed  from  the  z-planc  representation  to  a  u-plane 
representation*  In  this  new  representation  the  integration 
contour  coincides  with  the  real  axis  of  the  u-plane.  In 
the  third  section  the  u-plane  integrals  arc  evaluated  by 
the  common  method  of  steepest  descents.  The  resulting 
series  is  shown  to  be  of  limited  value.  In  the  next  section 
these  integrals  are  evaluated  by  a  modification  of  the 
method  of  steepest  descents  developed  by  Pauli.  The 
resulting  asymptotic  scries  is  compared  with  the  Pauli 
series  and  some  of  its  advantages  over  that  series  are 
pointed  out.  It  is  then  demonstrated  that  this  series  can 
be  rearranged  and  put  in  the  form  of  the  Oberhcttingcr 
scries*  In  the  final  section  the  equations  presented  in 
this  chapter  are  summarized. 

The  Saddle  Points  nud_  Poles  o_f  the  Integrand 

The  saddle  points  of  f(z)  arc  found  by  solving  for 
the  roots  of  the  equation 

f(z)  =  -j  ein  z  n  0  (P*0 

These  roots  arc  z  =  1<TT,  where  W  is  zero  or  a  positive  or 

c 
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negative  integer.  f(z),  therefore,  possesses  an  infinite 
number  of  saddle  points,  equally  spaced  alone  the  real 
Axio  of  the  z-planc.  In  view  of  the  locations  of  the 
caddie  points  and  the  behavior  of  exp  for  large 

values  of  Inz,  it  is  appropriate  to  create  a  closed  path 
of  integration  C,..,  consisting  of  c’  -  c,  SOp  ,  and  SDP^. 
This  path  is  shown  in  Figure  4.  SDP^anS  SDP^.  represent 
the  steepest  descent  paths  through  the  saddle  points  at 
t,  «  -*|fand  z  =TI“|  respectively. 

It  follows  from  the  Cauchy  residue  theorem  that 


*■  -it 


t  s  krf(z)  , 
(z)c  dz 


SDP 


F(z)ckrf,lZ)  dz  + 


IT 


(25) 

F(z)  has  poles  at  the  points 

t  e  -  ^  +  2nT/K  ,  (26) 

trhcrc  K  is  zero  or  a  positive  or  negative  integer.  These 
poles  lie  oft  the  real  axis  of  the  z-plane,  and  their 
positions  on  this  axis  are  determined  by  the  values  of  n,$*, 
and  ft 1  * 


■ 

2lTj2(Rc  siducs  of  poles  enclosed  by  Cy)| 


When  2n1i’IJ ) £  TT  ,  the  pole  is  enclosed  by 

®nd  its  residue  nust  be  included  in  the  evaluation  of 
Kkrjp,  n)  as  is  indicated  in  Equation  25.  This  residue 
can  be  found  by  first  writing  the  integrand  of  Equation  21 


- . . . . . . — 
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Z’PJLAWE 


rc  h*  Closed  integration  contour  in  the  corip  lo 
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tiD  the  ratio  p(z)/q(z),  v.hcrc 


jkr  cos 
cu 


z 


(27) 


and 


HM  =  '•ffjn  ein  (  — ) 


(28) 


Since  p(z)  and  q(z)  are  analytic,  R,  the  residue  of  their 
ratio  at  z  =  -  £3  +  2nini,  is 


n 


or 


p  (  7. ) 

q'(a) 


-  -e 


+  2mW 


(29) 


3__  jkr  coo  (-{3+  2TTnlO 
c  hij  c 

The  reddue  contribution  to  Gp  from  the  pole  at  z 

%  , 

io,  therefore, 


(2/0) 
-p+  2r/ifJl 


C*i,!r  cos  (“P+  2n‘,N)  u[if-|-p+  2nin;j]  , 


rherc  U 

* 

represents  the  unit  step  function. 

• 

f 

0  if  1  <  0 

«<t)  rJ 

1  if  t  >  0 

(31) 

l 

I  if  l  s  0 

The  last  requirement  in  the  above  definition  of  U(t)  accounts 
for  the  fact  that-  the  residue  of  the  po3c  has  the  Cauchy 
principal  value  v/hen  the  pole  lies  c>r»  the  integration  contour. 
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The  above  contribution  to  Gp  describes  a  plane  wave  propagation 
in  a  restricted  region  of  space.  If  N  =  0  and  (3  -♦-4  *  the  residue 
contribution  to  Gp  has  the  form 

eJkr  cos  <M<> 

This  term  describes  the  plane  wave  field  incident  on  the  wedge.  The 
angular  space  illuminated  by  this  wave  depends  on  the  angle  of  in¬ 
cidence^'.  The  illuminated  region  is  described  by  7 T~  -<p'l  >0. 

The  transition  between  the  illuminated  region  and  dark  region,  that  is, 
the  shadow  boundary,  is  defined  by  7T  -  |  ^  =  0.  If  c*.  -  7T  <  <71, 

or  ife*<7T»  no  shadow  boundary  can  occur.  When  a  field  point  falls  on 
the  shadow  boundary  the  pole  associated  with  the  incident  field  falls 
on  a  saddle  point.  If  -7T the  pole  coincides  with  the  saddle 

at  z  «  -7T,  while  if  0’>  7T  it  coincides  with  the  saddle  at  z  =  +7T. 

All  of  the  other  residue  contributions  to  Gp  describe  reflected 
plane  wave  components.  The  angular  regions  of  space  illuminated  by 
these  reflected  fields  are  described  by  7T  -  |2rtTTN  -  ($+^')|  >  0, 
and  the  boundaries  of  these  fields  by  7T  -  j2nTO  -  +  *  o,  except 

when  N  =  0  and  .<P  -  If.'.  Cl  early,  when  the  field  point 

falls  on  a  reflection  boundary  the  pole  associated  with  that  field 
component  coincides  with  a  saddle  point. 
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The  polec  of  F ( z )  play  two  roles  in  the  determination 

of  the  asymptotic  expression  for  G  .  First,  they  determine 

r 

the  plane  v/avc  components  of  the  field  through  the  residue 
contributions  discussed  in  the  previous  paragraphs. 

Second,  they  strongly  influence  the  form  of  the  asymptotic 
expression  for  the  integrals  in  Equation  25*  For  these 
N reasons,  it  is  important  to  have  a  complete  description 
of  the  pole  configurations  of  F(z).  In  particular,  the 
Index  N  of  the  pole  nearest  the  saddle  points  of  the 
integrals  in  Equation  25  must  be  known. 

Figures  5  and  6  illustrate  a  method  for  describing 
the  possible  pole  configurations  which  can  occur  as  ^  and $4 ' 
ore  varied;  o'  is  assumed  to  be  fixed  at  ICO  degreec  in 
this  example.  Figure  5  illustrates  the  possible  polo 
configurations  for  p  =  Figure  6  is  drawn  for 

pc  $5+  The  rectangles  denote  the  range  of  Rez  over 

which  the  polec  move  as  ,  the  angle  of  observation, 
varies.  The  right-hand  side  of  each  rectangle  coresponds 
to  e  0,  and  the  left  to  <p  =  oi .  Each  rectangle  corresponds 
to  a  fixed  value  of  N,  the  pole  index,  and  the  angle 
of  incidence. 

The  angular  regions  illuminated  by  the  incident 
and  reflected  fields  can  easily  be  determined  from  this 
diagram  as  can  the  locations  of  the  field  boundaries. 


/ 00 


\  rXfi,Ur0  ?'  Th°  allov,od  ranC*5  for  the  poles  of 
\M  for  various  values  of  ,*».  p  B  and  *  e  i00( 


j  C — >  /c?0 


r(y)  forU^*  The  allov/ed  ranges  for  the  poles  of 
1  J  f°r  Vari°US  Vftlueo  o^*.  /*  -  ind  oi  =100°. 


The  illuminated  regions  are  defined  by  the  valuer,  of 

• 

<f>  %  %  and  N  for  which  the  poles  lie  between  the  vertical 

lines  at  z  =  +  If.  The  field  boundaries  occur  for  the 

values  of  <*>  ,  f 5 ' ,  and  N  at  which  the  poles,  coincide  with 

* 

these  lines. 

It  is  important  to  observe  that  various  poles  can 
lie  near  the  c addle  points  at  z  =  +TTi  depending  on  the 

i 

value  of  $  and  ($',  In  Figure  6,  for  example,  if  ^5*  =  100 
degrees  and  =  100  degrees,  the  pole  N  =  0  falls  near 
the  saddle  at  z  c  -TTi  and  the  pole  N  =  2  lies  near 
t.  c  +TT*  On  the  other  hand,  if  s  0  and  $$  =  0,  the 

pole  H  «  -1  is  near  z  =  -  If while  the  pole  H  =  +1  is 
near  z  *=  +  TT* 

* 

Figures  7  and  8  picture  the  possible  pole  configura¬ 
tions  for  a  *i0  degree  wedge.  This  example  illustrates 
that  as  w  decreases  toward  zero  an  increasing  number  of 
poles  will  fall  within  the  region  -TV  £  Rez  £  TT,  and  that 
higher  order  poles,  that  is,  poles  described  by  larger 
Values  of  N,  will  lie  near  the  saddle  points. 

loiter  in  the  analysis  it  will  be  important  to 
recognize  that  the  pole  nearest  the  saddle  point  at 
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r,  =  ilT  is  described  by  the  N  which  most  nearly  satisfice 
the  equation 

2nlrN  -  (  ^  ±  tf**)  =  iir  •  (52) 

The  next  step  required  in  the  derivation  of  the 
asymptotic  expression  for  Gp  is  the  evaluation  of  the 
two  integrals  in  Equation  25*  that  is*  the  asymptotic 
evaluation  of 


and 
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7««nn» format. Ion  to  a  New  Complex  Plano 

In  order  to  evaluate  thecc  integrals,  it  is  moot 

convenient  to  transform  variables  to  a  complex  plane  in 

«,hich  SDP  coincides  with  the  real  axis.  This  new  complex 
••  II 

plane  will  be  called  the  u-plnnc.  The  transformation 
equation  relating  the  z-  and  u~plancs  is  found  in  the 
following  way.  For  values  of  z  on  SDP 

i'lf 


Hef(fc)  £  Kef(tfr) 


(55) 


Inf(z)  =  Inf(+TT)  . 


(56) 


Therefore,  on  SDP  ,  f(z)  can  be  written 

*Tf 

tM  =  f(ilr)  -  u2  , 


(5?) 


when  u  iG  real.  Substituting  Equation  22  into  27  yields 
the  transformation 


J  con  z  =  -  u4 


(58) 


4  TT 

•i  1  J* 

w  c  ie  v^cos  ?; 


(59) 


The  choice  of  ;:5.gn  must  now  be  made.  The  equation  for 


tho  steepest  descent  path  SDP^.  is 
too  x  cosh  y  =  -1 


(4«0) 


>6 


In  the  vicinity  of  z  =TF,  thio  reduces  to 

*  -tr  «  y  . 

•  • 

a  straight  line  making  an  angle  of  45°  with  the  x-axis. 
In  thic  vicinity  z  can,  therefore,  be  expressed  aa 


(41) 


%  e  TT  +  se 


if 


(<t2) 


where  c  is  the  dictancc  measured  along  the  path.  The  value 
of  c  io  zero  at  z  sTT,  positive  for  Imz  >  0,  and  negative 
for  Iraz  <  0.  It  follows  from  Equations  39  and  42  that  for 
%  in  the  vicinity  of 


<*» 

When  z  >  0,  u  must  be  positive  if  integration  along  SDP 
io  to  correspond  to  integration  from  -  oo  to  +oo  along  the 
roal  u-axic.  The  minus  sign  shou3d,  therefore,  be  used  in 
Equation  39«  The  transformation  which  maps  SDF^  onto  the 
real  axis  of  the  u-planc  is,  therefore, 

*  TT 

u  ts  VPcos  |  (W 

It  can  be  shown  in  a  similar  manner  that  this  transformation- 
fclco  caps  SDP^  onto  the  real  axis  of  the  u-plane* 


Figure  9  pictures  the  strip  of  the  s~plane  defined 


by  0  £  Rea  ^  2Tf  .  SDP.|r  is  also  shown  here.  Figure  10 

•  • 

ohows  how  the  various  regions  of  this  strip  map  onto  the 
u-plane .  Other  strips  of  the  z-planc  map  onto  other  Ricnann 
r.hccts  of  the  complex  u  representation.  The  strip  defined 
by  -211  <  Rcz  <  0  maps  onto  a  sheet  adjacent  to  the  one 
pictured  in  this  figure.  The  small,  lettered  circles  on 
these  figures  illustrate  the  manner  in  which  the  poles 

* 

of  F(z)  map  onto  the  u~plnnc.  Jiotice  that  the  poles  in 
the  u-plane  move  to  the  origin  .heir  counterparts  in 

the  z-planc  move  toward  the  saddle  *  z  =  TT. 

Applying  the  transformatir  defined  by  Equation  44 
to  Equations  33  and  3 yields 


Itl,(kr<  n) 


(*•5) 


where 


c(M)  =  Cot(<tt_»M) 

and 


dr, 

du 


r. 


TT 


.  z  ( u  ) 
cxn  -r~ 


(46) 
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Figure  9*  A  portion  of  the  strip  of  the  complex 
r.-planc  defined  by  0  ^  Rez,  £  2  TT.  The  dotted  line  denotes 
SDIVr  and  the  small  circles  represent  poles  of  F(z). 


U- PL  A  WE 


Figure  10.  The  portion  of  the  complex  u-planc 
c< rrcrpomJinjj  to  the  portion  of  the  s-plar.c  shown  in 
Figure  9«  Several  of  the  z-:0.nne  boundary  liner;  arc 
Indicated.  The  circled  nur.V  ro  indicate  corresponding 
ffCionr.  in  the  tv.o  planes.  coincides  with  the  real 


u— s. 
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Thuo  Equation  *<5  can  be  written 


I  (kr,  0.  n) 
ITT 


-J(kr  ♦  ^  (  cot 

t  \  2n 

1T2^n  I  -VI"»(u I 


-kru  .  , 
e  du 


(<i8) 


n) 


-J(kr  ♦  £) 

i 

#n’2/?n 


G(u)e‘**cru  du 


(*♦9) 


with 


0(u)  = 


sin 


04  z(u) 
2n 

z(u) 


(90) 


This  completes  the  transformation  of  I  (hr,  (3,  n) 

ITT  1 

from  the  z-plane  to  the  u-planc  representation.  In  the 
remaining  ccctionc  of  this  chapter,  this  intccrnl  will 
bo  evaluated  by  the  method  of  steepest  descents  and  the 
Pauli-Clennow  modification  of  thin  method. 


Evaluate  on  of  T^(kr,  0  ,  n)  by  th  e  Ke  thod  of  ft  tec  pert  Percent 
In  this  section  Equation  ^9  io  evaluated  by  the  ordinary 
cethod  of  steepest  descents.  In  this  evaluation  it  is 
orsumed  that  none  of  the  poles  of  F(u)  lie  rear  the  origin 
of  the  u-plane.  Thic  is  equivalent  to  aesuninc  that  none 
of  the  poles  of  F(z)  lie  near  z  =  4tT»  Physically  thin 


assumption  inplcs  that  the  diffracted  components  of  the 


field  are  being  evaluated  at  points  well  removed  from 
shadow  and  reflection  boundaries. 

If  G(u)  has  no  poles  near  u  =  0,  it  can  be  expanded 
in  n  Maclaurin  series  having  a  radius  of  convergence 
extending  from  the  origin  to  the  first  pole*  If  this 


expansion  is  made,  and  the  result  substituted  into  Equation 
1,9.  I  (hr.  (?>,  n)  takes  the  forr.i 


I  (kr 
tit 


-iOtr  +  If)  ^ 

•  n)~"  2ji  CU  (f’  n)  ' 

B  =  0 


where 


0 

tV) 


(52) 


The  positive  superscript  on  the  coefficients  of  the  Maclaurin 
series  indicates  that  the  derivative  of  G(u)  is  to  be 
evaluated  at  u  c  0  on  the  Ricnann  sheet  corresponding  to 
tho  si-plane  strip  0  £  Rea  £  2TT.  The  minus  sign  denotes 
that  the  evaluation  is  to  be  made  on  the  u -plane  Rieuann  . 


choct  corresponding  to  the  r.trip  -2TT  <£  Re 7.  <  0. 

The  odd  terns  of  Equation  integrate  to  zero 
co  it  can  be  written 


—  CO 


Integration  of  this  series  reduces  it  to  the  following 
form 


-j(kr  + 


-jwir  +  jjv  p 

n) - V572T— LC*  (f3’  "Vi?'  + 


*  £  4  «?■  »> 


1  •  3  •  5  •  •  •  ( ?r.i  -  1 ) 


Ut:l 


•  — ‘-ti 

»  ♦  2  I 

/» _ N  <•  J 


(5*0 


<kr)  c 

or,  in  more  compact  notation, 

-j(*r  ♦  jjo  .  f?,„  t  1, 

V:r,f''  n)~‘  VzTS'n  )  CL  <0-  n> - — T 

■nVo  <to>  2(55) 


hZ 


°p  »  I_K(kr<  ^  =  J*  -  f*' «  »)  ♦  I+1rO:r,  (3  =  f*\  n)  ± 

*  Wkr’  P=  **+S*'’  n)  + 

*  e*kr  C0G  * 

*  plane  wave  terms  corresponding  to 

reflected  field  components#  (59) 

•  #  * 

It  follows  from  the  derivation  presented  in  this  section 
that  co  long  as  the  field  point  at  which  is  being  evaluated 


doc/3  not  lie  on  or  near  a  shadow  or  reflection  boundary 
Gp  can  be  written  as  follows: 


jhr  cos  |„j_  j5,|) 


•f 


♦  plane  wave  terns  corrocpondinc  to  reflected  field 
components.  (60) 


This  equation  ic  identical  to  the  asymptotic  expression 
for  Gp  developed  by  Oberhettincer  for  the  case  in  which  the 
field  point  is  well-removed  fro n  shadow  and  reflection 
boundaries.  If  only  the  first  tern  ic  retained  in  each 
of  the  cur.nationc  in  this  equation  it  reduces  to  the 
familiar  fern  for  Gp  used  in  the  Geometrical  theory  of 
diffraction,  that  ic, 


44. 

GP  =  0  “  ^*)  1  D(n,  9^  ■»  ' ))  + 

♦  cjkr  cos  (Ji-Jf)  u(1)._  + 

+  plane  wave  terms  corresponding  to  reflected  field 
components,  (6l) 


where 


D(n,  i*  i  <?>') 


is  the  edge  diffr^ic tion  coefficient. 

It  was  pointed  out  in  the  examples  on  pages  29  through 
34  that  when  f*,  the  angular  coordinate  describing  the  field 
point,  approaches  a  shadow  or  reflection  boundary  a  pole 
approaches  at  least  one  of  the  caddies  nt  z  s  +11*.  In  such 
a  situation  the  corresponding  u-plane  pole  will  move  toward 
u  *=  0.  This  will  cause  the  radius  of  convergence  of  the 
Maclaurin  series  used  in  the  evaluation  of  Equation  49 
to  shrink  to  zero,  and  the  coefficients  of  that  series  to 
grow  without  bound.  This  renders  invalid  the  asymptotic 
series  expression  for  at  least  one  of  the  first  four  terms  of 
Equation  59#  A  method  must,  therefore,  be  developed  for 
evaluating  Equations  35  ar.d  54  for  the  case  in  which  n  pole 
of  F(z)  lies  near  the  saddle  point.  One  such  method  is 
discussed  in  the  next  section. 


Evaluation  of  I  (hr, ft,  p.)  by  the  Pauli-  Cl  err*  or  Modified 
Method  of  Steeliest  Descents 


In  thic  section  Equations  and  J>h  are  evaluated  for 
the  cane  in  which  a  single  pole  lien  near  z  =  -If  or  z  -  -i <TT> • 
This  evaluation  utilizes  the  basic  ideas  employed  by  Pauli 
but  yields  a  more  general  result.  The  approach  is  somewhat 
easier  to  follov/  than  Pauli's,  and  the  result  is  in  a  more 
convenient  and  useful  form  since  the  confluent  hypcrgeoi.ietric 
functions  are  not  introduced. 

The  essential  steps  in  thic  evaluation  arc  the  following 

a)  The  integrand  of  the  equations  is  first  written  as  the 
product  of  two  functions.  One  of  these  describes  the 
pole  near  the  saddle;  the  other  is  analytic  at  the 
origin  and  at  the  pole  in  question. 

b)  The  integral  is  then  transformed  to  the  u* plane 
representation.  The  analytic  portion  of  the  integrand 
is  expanded  in  a  Maclaurin  ncricc,  and  the  order  of 
summation  and  integration  is  interchanged. 

c)  The  resulting  series  ic  then  integrated  term  by  term. 

This  process  ic  described  in  detail  for  the  evaluation 
of  I  (krf  p,  n)  for  the  ease  in  which  the  pole  of  F(z)  at 
z  r.  -^is  near  the  caddie  point  at  z  e  -Tf.  The  results 
ore  then  presented  for  the  evaluation  of  n) 

fox'  a  pole  near  z  =  +TT.  The  technique  used  in  thic  section 


fails  for  snail  values  of  n. 
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It  ic  convenient  to  begin  our  diccuccion  with  Equation  33 

The  integrand  of  this  equation  hac  poles  at  z  =  -  p  4  2nTni  • 

•  • 

So  long  an  n  5.o  somewhat  greater  than  1  only  the  pole  at 
z  c  -  p  can  lie  near  the  naddlc  point  at  z  =  -IT.  A 
technique  for  evaluating  I_^(kr ,  ^ ,  n)  for  this  ease  will 
now  be  described*  If  the  integrand  of  Equation  33  ic 
multiplied  and  divided  by  (coo  z  -  coop),  then  it  can  be 
written 


•“  ( Hf)'1" 


:r  cos  z 


n(z) 


coo  z  -  cos 


c 


jkr  coo  z 


(62) 


where 

•  •  *  *  •  ' 

H(z)  C  cot  ^  cos  z  -  cocf)  •  (63) 


H(z)  is  analytic  at  and  near  z  =  -(*,  and  H(-f)  =  2n  cin^  • 
The  information  concerning  the  pole  at  z  s  «p  ic  contained 
now  in  the  function  (cos  z  -  cosp)*’^.  If  Wkp-P- n) 
is  now  transformed  to  the  u-planc  representation  it  takes 
the  form 


n) “  ■ 


-J(Kr  +  ]■•) 


OO 


j2/2'lTn 


.. /  %  -kru 

- d« 

4  ja 


(64) 


-00 


v/hcre 


H(u)  = 


(coo  z(u)  -  coo 


(6S) 


&  s  1  4  COS  p  (66) 

and 

P 

j(u'  4  ja)  =  cos  z(u)  -  cos  £  .  (6?) 


ll(u)  is  analytic  in  the  vicinity  of  the  origin  so  lone  as 
only  the  polo  at  7.  =  ~p  is  near  7.  =  -/Tf ,  and  can  be  expanded 
in  a  Maclaurin  series.  Therefore,  Equation  Ch  can  be 
written 


Z^Oir,  p,  nU 


-j(hr  4  ^) 


where 


s  i  a"  * 

,  n)  =  • — - l(u) 

tni  dum 


u  «  0 


(68) 


(69) 


A8 


All  of  the  odd  terms  of  this  scries  vanish.  The  remaining 
terns  can  be  expressed  as  combinations  of  well-known 
functions,  namely  the  Fresnel  integrals  and  inverse  powers 

of  kr* 

The  evaluation  of  this  series  is  no w  described.  The 
integrals  arc  treated  first,  then  the  coefficients  arc 
discussed.  The  integral  in  the  first  term  of  Equation  68  is 


icoo 


CO 

.  K 

f  u  <•  ja 

-CO 


du 


(70) 


where  X  =  kr.  This  can  be  evaluated  by  first  differentiating 
with  respect  toX»  integrating  on  u,  then  X  ,  and  finally 
•  again  on  u.  The  details  of  this,  procedure  arc  as  follows. 
Consider 


© 


-Mv  +  ja) 
(u*  ♦  ja) 


du 


-oo 


(71) 


Differentiation  with  respect  to  X  followed  by  integration 
on  u  yields 


CO 


(72) 


Integrating  now  on  X  results  in 


(75) 


The  first  integral  in  this  equation  can  be  evaluated  by 
integrating  over  the  closed  path  C  shown  in  Figure  11. 


Since  the  integral  over  the  semicircular  part  of  C  vanishes 
as  R  approaches  infinity, 


50 


CO 

Sdu 
u2  +  ja 


-  co 


<  7r  <  nr 

-j  r  ir 

.  (u  +  C  /f?)(u  -  c  .  /o') 


2Tf j  ^ KcciducG  of  poles  enclosed  by  C*|* 


The  residue  of  the  pole  at 
,  Tr 

-i  JT_ 

u  *s  -  c  Va^ 


<-«/!> /’  V‘ 


and,  therefore, 


rm 

It111-' 

\  u  ♦  oa 
-oo 


/T 


if 


Substituting  this  recult  into  Equation  73  yields 


(7'0 


(75) 


I  (kr)  . 

o 


ir1 


„  jOa-n  ♦ 

V-e 


(76) 


The  integral  in  this  equation  can  be  expressed  as  a  complex 

TT  2 

Fresnel  integral  by  means  of  the  substitution  n^’  r.  ~  t  • 
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The  result  is 


I  (kr) 

o 


Try^f3kra 

/«' 


where 


(77) 


Thic  can  oloo  be  written 


(78) 


The  evaluation  of  the  hichcr  order  integrals  in 
Equation  68  is  noy/  quite  straightforward*  The  integrand 
of  the  next  even  tern  can  be  expanded  ar.  follows: 


u 


'  2 

u*  ♦  ja 


1  - 


2  4 

u  ♦  jo 


(79) 


Therefore, 
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( 


hr  / 


jal^vhr)  • 


(80) 


Tho  integrand  of  the  tern  corresponding  to  n  *  h  can  bo 
written 


n*  2  ,  (-In)2 

•5 -  s  u  -  ja  4  "y* - 

u  ♦  ja  u  4  ja 


(81) 


eo 

.00 


/iT 


2 

u  4  ja 


•  00 


2(ki-) 


1 

2 


ja  4  (-ja)2  I0(hr) 

<kr)2 


(82) 


Tho  integral  of  the  general  2m**1  term  can  be  written 


2ra  .  2 

u  _  -kru  . 

—  e  du 


u2  +  ja 


•eo 


■ 

ZX.3,5>,..(p(n  -  );)  ..  p/iT  # 
2<»  *  k><kr>"  “  k  +  * 


k  ■  1 


(-Ja)k  '  \  (- jn)n  Io(kr) 


(8 » 


or 1  in  more  compact  form, 


5) 


$4 


It  foil  owe*  from  Equations  65,  6y,  69 •  and  50  that 

(86) 

Thio  yields  the  following  rolations: 

B©  (pt  n)  e  -aC©  (^,  n)  (87) 

and 

B-o  (f,  n)  .  JC *(n  _  1}  (p,  n)  -  f,  „)  .  (88) 

Since  1J(u)  is  analytic  when  z  =  -p,  all  of  these 
coefficients  are  finite  so  long  as  only  the  pole  at 
*  e  -  p  approaches  the  saddle  at  z  r  -IT.  In  particular , 
they  arc  finite  when  p  fi'TT.  The  follov/ing  analysis  reveals 
tho  value  of  these  coefficients  for  this  value  ofp. 

H(u)  can  be  written  ac  the  cun  of  an  even  and  odd  function 
of  u,  that  is, 

H(u)  .  11® (u)  ♦  ll°(u)  ,  (89) 

where 


1  JJ™ 
(2m)i 


{.IOCuXu2 


H*(u)  c  PfM)  *  "H 

end 


(90) 


H(it)  -  H(»u) 
2 
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n  (u) 


(91) 


In  order  to  find  the  form  of  lr(u)  and  It(u)  it  ir.  nost 
convenient  to  r.ahe  the  change  of  variable 

»(u)  ciO(u)  -IT 


(92) 


oo  that 


*  ,tr  ir 

II  m  -C  V?  COO  £  e  -C  Gin 


(95) 


The  nllov/cd  values  of  U)  arc  restricted  to  the  chndcd  ctrip 
of  thctO-plnne  chcv:r,  in  /inure  12.  Therefore,  a  chance 
in  the  cicn  of  to  corresponds  to  a  change  in  the  oicn  of  u. 
This  permite  Hc(u)  and  H°(u)  to  be  written  ac  follows: 


ir(u) 


conO(u)  <  cor- ft 

^  .  tofu ) 

2  coo  — 


4  COt 


-rt(u) 


and 


U°(u)  = 


«{u)  «  p  -ir 
cot  - X - ♦ 


i**] 


:oo*d(u)  ♦  cos  ft  r  .  *^(u)  «  ft  -tT 

T - C3HD -  cot - S - 

2  coo  -g-  |_ 


(9'0 


-vX'ii)  *  (S  -TT  *( 


(95) 


cot 


How  when  p*-*  ♦TT  «  H°(u)  e  0  and,  thereforo,  all  of  tho 
even  coefficients  of  the  Ihclaurin  ccrico  expansion  of 


H(u)  vanish,  that  its, 
B;n  (IT,  n)  *  0 


(96) 


Figure  1 ?.•  The  coupler  O- plane#  The  values  of  u> 
allowed  in  Equation  93  arc  restricted  to  the  shaded  regions 
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For  thio  caoc  of  p  «  IF  Equation  85  reduces  to 


-J(kr  ♦  ]j0 


W'r.'ff.  n)~  "jg/p,,  tr  •  no(1r'  n)  Io<1:r> 


(97) 


Thio  firct  terra  Bust  be  retained  since  the  pole  in  I  (kr) 

o 

•t  ?  ntT  juot  cancels  the  aero  in  Bo(p,  n)  to  yield  a 
finite  value.  That  value  in  found  by  writing  (3  «  TTi  € 
and  taking  the  Unit  ao  C-*-0.  The  details  arc  shown  here. 


i€ 


B0(ltie*  n)  C  -  (1  ♦  confir  1  £*]  )  cot  jj- 


X  7  r.G. 


(98) 


»  ,  ir /?  -if 

I0(kr) *  —  e 


•o 


(99) 


11m  I  ^(kr,  IT  t  6  1  r.)  *  t  v-  e“^r  (100) 

6 -VO  2 

The  plus  ticn  applies  when  p>TT,  and  the  minus  ci£n  when 
p<TT •  Thio  discontinuity  in  the  diffracted  field  juct 
compensates  for  the  discontinuity  at  the  boundaries  of 
the  incident  and  reflected  fields,  producing  a  continuous 


total  field  there. 


The  value  of  Equation  100  directly 


on  the  field  boundary  is  zero.  The  Cauchy  principal  value 

of  the  associated  reniduc  contribution  in  £  «  no  k&o  ohor/n 

•  • 

in  Equation  31. 

Th©  pole  nt  z  *  -  p  v  -  (fi-  $<•)  con  lie  near  the 
toddle  point  nt  z  ■  4 IT •  V/hon  thin  situation  occur© 


aunt  be  evaluated  in  the  manner  jur.t  illustrated.  This 
•valuation  yields 


p,  n)~ 


•J(V,r  4 
J 


co 


m  *  0 


Z[(n  -  k  ♦  IX-ia)*  ’  1 

(u  -  k  ♦  h) 

(kr)  2 


♦  (-ja)n  I0<kr) 


With 

lira  I  (kr,  ir  i£  ,  n)  a  *  | 
C+O  2 


(101) 


(102) 
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#*<p,  n>  .  -aC*((>,  n)  ,  <10}) 

and 

B^(P,  n)  -  JC;(n  .  ^(p,  n>  -  .c;o(p,  n)  .  (10*,) 

Figures  5*  7 *  and  8  illustrate  that  the  poles  of 

F(s)  located  at  z  ■  -jj*  ?Y> ftn  can  lie  near  the  caddie  point 
at  a  •  -TT  when  II  is  a  negative  integer*  When  this  situation 
occurs  I^(kr,  p,  n)  con  be  evaluated  by  the  Method  used 
for  the  case  in  r/hieh  the  pole  at  z  «  -p approaches  that 

oaddla*  However,  the  method  mist  be  nodified  slightly  to 

« 

account  for  the  fact  that  tho  pole  is  described  by  arbitrary 
negative  N*  This  modification  entails  simply  replacing  the 
term  (cos  z  -  coop)  of  Equations  62  and  63  by  (cos  z  • 
cos  (-p*  PritfK))*  This  substitution  Makes  ii(»)  analytic 
at  *  *  -  p  ♦  2ritni  and  allows  expansion  of  J!(u)  in  a  Eaclaurin 
series*  This  substitution  puts  Equation  66  into  the  form 
a  «  1  ♦  coo  (- (3  ♦  2riTfII)*  All  other  cqoatior.3  of  the 
analysis  retrain  unchanged  in  forra*  Therefore,  when  Gp 
lo  to  be  determined  at  field  points  for  which  -  (1  ♦  2rifl‘h* 
is  near  -IT,  Equation  8 rj  should  be  used  to  evaluate 
Wk*-P-  »>•  In  Equation  83  and  the  allied  equations, 
that  is,  87  and  88 ,  a  s  1  ♦  coo  (-  p  ♦  2nTT,0. 

It  also  can  be  shewn  that  when  Gp  is  to  be  determined 
at  field  points  for  which  -  p  +  2r.VN  is  ncarlT,  (sec  Figure  fi, 
for  example)  Equation  101  should  bo  used  to  evaluate 


I^Ocr,  p,  n).  In  thin  equation  a  «  1  +  coo  (-  <  2nTfJ.'). 

In  this  cane  II  in  a  positive  integer. 

If  n  io  clone  to  */cro,  that  in,  if  the  exterior  angle 
between  the  facen  of  the  scattering  wedge  io  a wall,  a 
lnrgc  number  of  polcc  will  lie  near  the  saddles  at 
*  «  iff*  In  thin  cane  the  rad 'mu  of  convergence  of  the 
MaclnurJn  series  for  l!(u)  wil  be  very  small  ant*  Equations 
and  101  will  fail* 

A  lumber  of  examples  ill  not  rat  lug  the  uto  of  the 
equations  developed  in  thin  ncction  are  diccuorcJ  in 
Chapter  IV.  However,  before  proceed ing  with  thin  the 
relation  between  the  expressions  developed  here  ard  thone 

given  by  Pauli  and  Obcrhettingcr  will  be  examined  further. 

• 

Relation  to  PmiH  *t  Arvrwtct*  c  Series 

Pftuli  wan  the  first  to  introduce  the  method  of 
•tccpcot  dccccnts  used  in  thin  work.  He  applied  this 
■ethod  to  the  evaluation  of  Gy  in  a  tanner  such  like  that 
Used  in  thin  study.  Pauli,  however,  .cr  forced  one  operation 
at  the  beginning  of  his  evaluation  which  caused  hin  resulting 
•aynptotic  expression  to  be  of  lees  general  applicability 
than  the  one  derived  here.  The  starting  point  for  I’auli'o 
derivation  v;ac  Sorter  fold  *s  nrr  a  tr.  which  io  given  by 
Equ.ation  IV*  Pauli’s  path  of  integration  for  thin  equation 


wao  closed  in  the  manner  oho v/n  in  Figure  4,  And  the  pole 
contributions  to  Gp  were  evaluated*  Thin  reduced  Gp 


to  tho  form 


°p" 


n,  ^  -  j<*)d*  ♦  ^  F(»,  n,  0-j$»)d»  t 


SDP 


SDP 


i  ^  F(s,  n,  p  ♦  ji * )d*/-  t  ^  FU,  n,  ^4  •  >d: 

SDP  _  SDP _ 


4  Residue  contributions  , 


(105) 


whore 


r(*.  n,  ?  t  s»’>  .  - 


1  r 


j  - 
*  n 


—5-,-  c*fcr  eoc  *  (106) 


By  anfcing  the  cubstitution  r.  *  00  -  If  in  tho  firct  ond 
third  integrelo  and  z  *  to «  ir  in  the  second  and  fourth, 
Psuli  reduced  Gp  to  the  sum  of  two  integrals* 


°P  *  21Qn 


SDP. 


coo  -*  •  COO 
V) 


e*Jl:r  coou> 

—  u>;t<»  rgrr dw* 


r&“""  1  zrfr 


•JV.r  cooto 

to «  Co 


r  du)» 


GDI* 


4  Residue  contributionc  • 


(107) 


4 


* 
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%  e^T  (01  -U  )  but  none  lies  near  z  *  -'tT(orTf'),  then  a 
single  pole  will  lie  near  u3r  0  and  Pauli's  method  of 
steepest  descents  yields  a  useful  result  for  Gp.  If, 
however,  one  pole  lies  near  z  r  ♦  IT  and  another  near  z  s  -IT, 
two  poles  will  lie  near  tO=  0.  In  this  case  Pauli's  method 
fails  and  his  asymptotic  expansion  for  Gp  in  not  valid. 

Pauli  reduced  the  generality  of  has  result  when  he 
combined  the  integrations  over  SDP^  and  SDP  into  a 
oingle  equation.  The  generalized  form  of  Pauli's  expres¬ 
sion  developed  in  thin  work  can  be  applied  to  cases  in 
which  polos  simultaneously  lie  near  z  =  +TT  and  z  =  -IT, 
since  the  integrations  over  SDP  and  SDP^.were  evaluated 
individually. 

Pauli's  asymptotic  expansion  for  the  integrals  in 
Equation  10?  is 


3  .  *TT  ( 

2tfjn  Cin  n  \ 


-  jkr  cog  to 


ir  io+  (>5 

cos  -  -  coo  - - — 

n  n 


SDP 


PB  = 


1  'T-Kkr-f) 
-zzz-sz-  sin  -  c 
n  n 


z 

n  b  0 


£(m  ♦  h  , 

j  • — r =~-  Ad  (p).S  (kra)  — — 

yir  n  (kl.)n 


(103) 


where 


S^kra)  c  (m  -  -^(kra)^  F^d,  -n  +  |  ,  jkra)  , 


(109) 


~Sm  2  .(?n)r  - 

c  d  cor.  to  +  cor(J 


(110) 


A  (a)  „  £ _ J2 _ _ c OT AX Mi 

A2.nT;  •  (2b7I  ,.,(2n)  'Tr  to  1 Q  to  UJ 

a  to  I  cog  -  -  cog  — “ —  cor.  £j  l 

10  ~  0 

and  a  =  1  +  coo  =  1  +  cor  (  +  *)•  The  functions 

F^d*  -rn  ♦  ^  ,  ^ikra)  are  confluent  hyper  Geometric  functions 
These  can  be  v/ritten  in  the  form 


7  1  1 

F^d*  -m  4  £  ,  jkra)  =  (n  -  -)(kra) 


2e;jkra 


hr  a 


dt  . 


(Ill) 


If  (as  Ott1  bar.  shewn)  the  integrals  in  this  equation  arc 
repeatedly  integrated  by  parte  Equation  108  can  be  written 
in  the  form  shown  below. 


.  cin  ?  -  j  (hr 

I*  JL  _  I* 

WB  n 


Vf>JB 


(-ja)k 


^  f(n  -  k  +  |.) 


(kr) 


n  -  h  +  ^ 


m  =  0 


k  =  1 


♦  (-ja)”  Io(kr)  . 


(112) 


Pauli's  expression  for  Ug  can  be  derived  by  combining 
Equations  85  and  101  for  a  common  value  of  N,  such  as  h  s  0. 
Combining  these  equations  results  in 

Up  **  I4^.(krf  n,p)  ♦  I-tT(kr%  n,p)  c 


-jO:r  -  jj) 
c _ 

7T/Fn 


(bL(?'  n)  +  BL(P» n)) 


B  c  0 


Z(-jtt)k  ■  1  [<n  -  k  ♦  |) 

(n  -  k  ♦  i) 

,  (kr)  2 


♦  <-ja)n  I0(kr) 


.  (1X3) 


It  can  easily  be  shorn  that 

Bk<?*  ">  +  Bin(P- 

- Z - i  oin  £T  A2n(P'  n)  5 

hence  e  and  the  relationship  of  the  Pauli  series 
to  the  generalized  series  is  established. 


Relation  to  Cborhottinrcr 's  Any not  otic  Series 

Oberhettinger  has  studied  the  diffraction  of  plane 
waves  by  wedges  which  have  exterior  angles  greater  than 
180  degrees.  He  showed  that  the  total  field  can  be 


expressed  as  the  sura  of  four  integrals  plus  a  number  of 
plane  wave  terns.  The  integrals  describe  the  diffracted 
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field  find  the  plunc  v/nvc  terns  describe*  the  incident  and 
reflected  fields;  as  in  the  solutions  previously  described, 
Oberhet linger *s  expression  for  Gp  is 

Op  n  I(kr,  ('IT-  (  0-  n)  +  I(kr,  (TT  +  (0~  0'))»  n)  ± 

i  I(kr,  (IT-  ($$+£•)),  n)  +  I(kr,  (7T  +  (0+0'))i  n)  + 

•  c;,kr  coa  (<4‘  J40  u(ir-|pi- fS-| )  ±  c3kr  coe  ({4*  ^  . 

•  OCTT-|yi+  s**l  >  ♦  eJkr  cos  (^i+>5,  -  2n1r)  • 

•  v(ir-|  ft*?  '  -  2ntr|)  (ii4) 


The  relation  between  this  integral  and  the  integrals  given 
in  Equations  33  and  3^  will  now  be  established.  Equations  33 
and  Jfh  have  the  form  . 


Replacing  a  by  io-1T  in  the  integral  over  SDP_^and  by 
in  the  integral  over  SDP^ yields 


t 


6? 


COGUi 


du>  (117) 


The  path  SDPq  Is  shown  in  Figure  13.  The  integrand  of 
Equation  117  can  be  written  as  the  sum  of  an  even  and  odd 
function  of  l*>.  The  integral  of  the  odd  part  over  SD?o 
vanishes,  leaving 


If  SDPq  is  deformed  so  as  to  be  coincident  with  the 
In tO-axis,  and  then  is  rotated  by  90  degrees,  by  replacing 
CO  by  jx,  Equation  118  takes  the  form 
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Reversing  the  order  of  the  argument  to  'ft  +  (3  and 
replacing  this  by  S  reduces  the  integral  to  that  usjed  by 
Oberhettinger,  that  is,  Equation  115.  The  following 
notntional  correspondence  exists  then  between  the  integrals 
used  in  this  work  and  those  used  by  Oberhettinger . 


I  (kr,  (pS-^'),  n)  =  I(kr,  T  -  (f4-j50,  n>  (120) 

I^(kr,  (&-  fi<),  n)  =  I(kr,  V  +  (<*-{*<),  n)  (121) 

+Tl 

I^Ckr,  (0  +  f*'),  n)  =  I(kr,  'if-  (  rA  +  ^  * ) ,  n)  (122) 

I+^(kr,  (f$+  f*1),  n)  s  l'(kr,  +  (  S*  +  n)  (123) 


‘  Oberhettinger  used  Watson’s  lemma  to  derive  an  asymp¬ 
totic  expression  for  Equation  115  valid  for  cases  in  which  , 
the  angular  coordinate  of  the  field  point,  does  not  fall 
near  a  field  boundary.  That  expression  is  given  in  Equation 

124. 
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with 


A*(S,  n) 


2Tl/2n 


cot 


'  2n  ' 


(125) 


A?($,  n)  =  -  —^7=  cot 

1  2fTn/? 


( * 

f  1  1  2  | 

r  $  \1 

V'2n  ' 

l*  1  i? coc  1 

^  2n  )  [ 

(126) 


and 


A®(5,  n)  c°t(5|)(|*(^-^7r)cSc2(|r)t 


1  /» 
♦  ^  CBC 

nH 


(4)1 


(127) 


This  coefficient  differs  from  the  form  Given  by  Obcrhc ttinger 
by  a  factor  of  one-half.  ObcrhettinGcr 's  coefficient  is  in 
error.  These  coefficients  have  the  same  form  as  the 
C*<fc,  n)'s  derived  in  this  work;  in  fact, 


*:«•  n)  =  Wr7?C°1?'  n) 

Ai(5,  n)  =  Fife  C2(P’ n> 


(1?8) 


(129) 


A2(5,  n>  °  gll/k  Vf*  n) 


(150) 
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and  in  general 


A°(£,  n)  = 

ID 


/  i\«  ± 

2TT/2h  C2ra(? 


(151) 


Equation  124  is,  of  course,  identical  to  Equation  55*  They 
were  derived  from  different  forms  of  the  sane  integral. 
Equation  55  was  derived  by  a  direct  application  of  the  method 
of  steepest  descents;  Equation  1 24  w as  derived  by  Watson's 
lemma.  Both  equations  fail  on  and  near  field  boundaries. 

In  order  to  derive  a  more  general  asymptotic  form 
for  I(S,  n),  valid  near  field  boundaries,  Obcrhct linger 
split  Equation  115  into  two  parts.  The  first  parr  was 
associated  with  the  diffraction  by  a  half-plane;  the  second 
described  the  modification  of  this  half-plane  diffraction 
term  caused  by  the  finite  wedge  angle.  The  first  part 
reduces  to  the  Fresnel  integral  and  an  infinite  series  in 
inverse  power  of  kr.  Oberhcttinger  applied  Watson's  lemma 
to  the  second  integral  to  form  an  asymptotic  series  in 
inverse  powers  of  kr.  Obcrhc ttingcr 's  series  is  shown 
in  Equation  132. 
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oo 


TT 


I($»  n) 


/v*  — 


e 


-jkr  co c$  +  j  £ 


/TT 


ogn 


S  \  c“jr  ar-  (122) 


(2kr)^  sin  |> 


a> 


.  e~jKr  '  3  4  ^  (J)-I(-  ♦  |)  n)  - 

a  s  0 

-*.<»]  — »-T 

U  m  +  tr 


(kr) 


where  X-(<S)  =  A  (£,  2)  +  A  (2V-«S»  2)  = 


in 


m 


m 


(-1)” 


2lT2  *  (sin  |)2n  +  1 


Although  it  is  not  mentioned  explicitly  by  Obcrhcttinger , 
«l«l>  then  S  raust  be  replaced  by  S  -  2<<in  Equation  132, 
If  S  ie  replaced  by  If  +  p  this  soriec  takes  the  form 


7? 


7} 


jkr  coop  ♦  J  f  /  2 

X(p*  n)/w  -  - - yjp - ccn  (inp)  V  dT- 

(krn)^ 


-Jkr  - 


iT\  5ni>  r^i 

/  ■  ♦  i  <!Tt/Jn 

La  (kr>  2  u 


■  ■  0 


p  1 

I coa  ?. 

T/5V’  4  ». 


<1»> 


the  correspondence  between  Cbcrhettinjer *n  expansion 
•nd  the  generalised  Pauli  expansion  will  now  be  investigated, 
and  it  rill  be  shown  that  they  are  cluply  different 
arrangcacnto  of  the  cone  series.  Consider  again  Equations  $5 
and  101.  Thcoe  equations  are  written  below  in  a  slightly 
Modified  forn. 


Au  before,  a  =  1  +  cog  (- p  +  PritfiO*  but  here  the  Fresnel 
integral  has  been  factored  out  of  the  summation.  The 
ounnaticn  multiplying  this  integral  can  be  written  in 
closed  fern  by  setting 


<-ja)n  ,•  (/TJ7?) 


-.vO  % 


eo  that  the  summation  has  the  form 


<ri 

o  =  0 


-'v  ?n 


Tliic  is  just  equal  to  the  even  part  of  J!(u)  (see  Equations  Ctjl 
89,  and  9;0  evaluated  at  u  =  s/ - ,  This  value  of  u 
corrccpondc  to  z,  *  -^+  2nTVN  or  to  u>:.-  ~p  +*\r+  2nTi1\ 

(ccc  Equation  92).  The  even  part  of  H(u)  can  be  therefore 
evaluated  as  follows: 


n*(u 


(cor;iO(u)  +  cor>(-fi  +  niffiO)  # 

~  _  to(uT" 
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coc - -- - 


ein 
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ui  ~  TT-  p  f 
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(w> 


Tj 

The  eecond  tern  of  this  equation  vanishes.  Application 
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The  factor  sgn  (TT±  ((3-  2riIfiO)  cnriscG  here  because 


The  Minus  sign  applies  if  -  2ninij>'n*j  otherwise  the 
plus  sign  applies. 

The  second  tern:  in  Equation  1^6,  that  is,  the  double 
©urination,  can  be  rearranged,  and  written  as  a  series  of 
coefficients  multiplying  terns  of  the  form 

i  <  • 

-M 

<kr)2  . 


That  is, 


m  «  1  .  k  =  1 


m  =  2 


77 


The  typical  coefficient  Sp(^,  n)  can  be  expressed  as  a 
finite  sum  by  making  use  of  the  series  form  of  the  even 

part  of  K(u)  as  shown  here. 

oo 


Sp^,  n) 


<-j=rPtl  2^  B*ra(p,  n)(-j.,)B  = 
m  =  P  +  1 

(-jar01  +  «/=3?  >  - 


n) (- ja) 


(1^2) 


The  final  sum  in  this  expression  can  be  simplified  using 
Equations  87,  88,  102,  and  10*4.  Substituting  these 
equations  into  the  sum  yields, 
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r 

£  4<e-  n)(-^)m  - 


m  n  0 


-<‘C  (aC  •*  j1"1  C«)  +  (-ja  C,,  +  a^C.  )  + 


7.  i  L  t  p 

+  (-a^C,  -  ja  CA)  +  ••••  +  (-j)  • 


P  +  1  ±  /  .NP  P  +  1  ±  yv 

a  ^2P  =  '“0^  a  ^2P  1 


eincc  nil  but  the  last  term  cancel .  Substitution  cf  this 
result  and  the  expression  for  the  oven  part  of  l!(u)  into 
Equation  lh2  produces 

Sp(^,  n)  =  <-jc)-(P  +  15  (-2n  coo  (f  -  nTT;i ))  ,  jC^,((3,  nXlW 

This  equation  can  now  be  inserted  into  Equation  l4l  and 
that  alone  with  1^0  and  I36.  These  substitutions  yield 


X+  (kr,p,  n)  ^  -  sgn  (TT+  (p-  2nTflO) 


jkr  cos  (p-  2nTT;;)  +  j  £ 
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(X'i5) 
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This  i g  just  Obcrhcttinger ’s  series  as  Given  in  Equation  125* 
but  here  written  for  arbitrary  wedge  angle.  Thus,  the 
Obcrhc ttinger  series  is  simply  the  Generalized  Pauli 
series  written  in  a  different  form. 

Summary  of  Equations 

In  this  chapter,  five  scries  expressions  for  Gp  have 
been  presented.  The  first  of  these  was  the  ci£enfunc tion 
series,  and  the  other  four  were  asymptotic  series.  These 
five  expressions  are  summarized  below,  and  their  range  of 
application  is  discussed. 


'Eigenfunction  Series 


j  ~  IT 

=  It  )  €«Ja(kr)c  n  2[C0S  n 
m*  0  n 


’ )  ±  cos  • 


G  <*♦*.)] 


(1V6) 


This  equation  is  perfectly  General  and  can  be  employed  to 
compute  Gp  for  arbitrary  wedec  angles  (*<) «  angles  of 
incidence  (f!1 ) ,  and  observation  coordinates  (r,  ff). 

The  four  asymptotic  expressions  for  Gp  all  have  the 
following  basic  form: 


u(ir~ 


IjM’I) 


Gp  “  c 


jhr  cos  ( $  -  ft 1 ) 


■» 


8o 

+  (plane  wave  terns  corresponding  to  reflected  field 
composes  to) 

*  ”*  n*  i 

t  n,  (^♦^’))  t  I^Oir,  n,  (f^+f*1))  (l1*?) 

The  first  tern  in  this  expression  describes  the  incident 
plane  v/ave  field.  The  second  tern  contains  the  plane  wave 
field  r#  rising  fron  the  reflection  of  the  incident  field 
by  the  faces  of  tlie  wedge.  These  fields  are  described  by 
terns  of  the  form 

bJkr  co°  iitm  -  f,))  udr -  1 2nu;:  -  (jJt  t<)\ ) 

The  value  of  K  and  the  range  of  existence  of  these 
higher  order  components  arc  most  easily  determined  from 
the  a- plane  pole  diagrams.  The  unit  step  functions 
appearing  in  this  equation  describe  the  locations 
of  the  plane  wave  field  boundaries.  Cn  the  field 
boundaries,  the  unit  step  has  the  value  cne-half. 

The  four  final  terras  in  Equation  1*1?  describe  the 
diffracted  component  of  the  total  field,  and  it  is  these 
terras  which  are  expressed  as  asymptotic  series .  Each  field 

t 

boundary  has  associated  with  it  one  of  these  four  terms. 


The  four  asymptotic  forms  for  these  terras  are  listed  below. 
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Unconnected  Abymptotic  Series 


-J(l:r  +  jj-) 


Itl(kr.  (*t  ji'),  n)  -  -  z/rnV  ■ 

f1  i  [?n  ♦  |) 

L-X  (kr )"  +  2 


o  =  0 


/  Tfl  (3  \ 
C«  c  cot  (  -3T-  ) 


(l'<8) 


(149) 


(150) 


•  (151) 


This  expression  con  be  used  to  conputc  the  four  final  terns 
In  Equation  1^7  so  lon£  as  the  valuer,  of  <f>y  <J> ' ,  and  n  are 
such  that  no  polos  lie  near  the  saddle  points  associated 
t/ith  the  evaluation  of  these  terns.  Tf  the  values  of 
s*.  $<\  and  n  arc  such  that  a  pole  lies  near  the  saddle 
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associated  with  a  particular  I  ,  Kouation  1^8  is  not 

tV  • 

valid . 

t 

When  hr  is  increased,  the  effective  distance  between 
the  poles  and  saddle  points  is  increased.  Therefore,  for 
given  values  of  f$,  <£'  and  n,  the  validity  of  Equation  l^S 
increases  with  increasing  kr.  Stated  in  physical  terns, 
Equation  1>«8  can  be  used  to  compute  the  diffracted  components 
of  Gp  at  field  points  at  which  the  product  of  the  angular 
distance  from  field  boundaries  and  the  radian  distance  from 
the  ertex  of  the  wedge  is  large. 


Generalised  Pauli  Asymptotic  Sorici 
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BIm(^±  n)  =  JC”(n  _  x)  ((jjSijd'),  n)  -ac|o  . 

( ( ±  * )  ♦  n) 


(15*0 


a  e  1  +  cos  ( 2n7fN  -  (^  +  ^*)) 


(133) 


I  (kr)  =  §f>ra  c‘jt  dt 

O  /TT*  \ 


(156) 


(kra)' 


N  denotes  the  s-planc  pole  nearest  the  saddle  point  used  in 
the  evaluation  of  I+  .  Equation  132  can  be  used  to  evaluate 
the  last  four  terns  in  Equation  1*»7  for  arbitrary  values  of 
$%  <f>x  and<<,  and,  in  particular,  for  values  of  ^  near  the 
boundaries  of  the  plane  wave  fields.  The  value  of  N  to  be 
used  in  each  tern  is  most  easily  determined  from  a  pole 
diagram.  The  accuracy  of  this  expression  increases  with 
increasing  kr. 


Pauli  Asymptotic  Scries 


1  (krtp»  n)  +  1  (hr,  ft,  n) 

-If  '  \ 


+TT 


IT  ...  'Til 
1  "°(  71  "  V 

A  'ir  /^‘  ” " c 


+ 

The  expressions  for  K-:|((3,  n)  arc  given  by  Equations  1DJ> 
and  15;i,  and  Io(l;r)  is  given  by  3^6. 

Equation  157  nay  be  used  to  compute  Ihc  final  terns 
of  Equation  1*»7  in  cases  in  which  a  pole  lies  near  cither 
the  caddie  point  at  z  rr'TT  or  at  z  =  but  not  near  both. 
Thic  equation  is  a  cun  of  two  series  of  the  form  given 
by  Equation  1^2  for  equal  values  of  and  II  -  0. 
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Generalised  Oberhet tinker  Sorico 


I  (kr,  p,  n)^  -  ocn  (IT  i  p  -  2nTi:).  • 


$kr  coo  (p-  Zrffll)  ♦  J  ^ 

1  7^ 


.-it2 


(kra) 


1 

2 


-J(kr 


•  o 


*  T>  V  Dp  *  ?>  f cgp<c.  «■)  _ 

/  P  ♦  i  L  n  T 

Ljl  (kr)  2 


P  •  0 


$  COS  (5  - 


IT  /S\ 


-  nTx) 

V  «  1 


05S) 


a  «  1  4  coo  (p-  2n1Tli) 


The  equation  can  be  used  to  evaluate  the  loot  four  tcrr.r» 
in  Equation  1*»7  for  arbitrary  vnluec  of  f**  andK  , 

The  accuracy  increases  rith  increasing  kr.  Thic  series 
io  ciraply  a  rcorrar.f.crent  of  the  ccnc**alined  Pauli  scries, 


CHAPTPK  1  V 


’;>;a;;i  iro 

Sorie  typical  numeric;  3  example;  arc  presented  here 
to  illustrate  the  tyre  of  results  which  arc  obtainable 
with  the  equations  (level  opt  d  in  the  previous  chapter. 

Some  conclusions  regarding  the  relative  merit  c»l*  the  ordinary 
Pauli  series,  the  f eaeral ixed  Pauli,  aeries,  and  the 
Oberhe ttinger  aeries  are  sur.i.iari sod  at  the  end  of  this 
chapter. 

In  each  example,  the  wedge  angle,  the  angle  of  incidence 
of  the  i>lane  wave  field,  the  ancle  of  observation,  and  the 
boundary  conditions  v.cro  fixed,  and  the  total  field  was 
calculated  ns  a  function  of  kr.  In  nil  the  examples,  the 
calculations  were  made  on  or  near  a  field  boundary  for 
relatively  small  values  of  hr;  in  most  cares,  0<  hr  <  lrj. 

It  is  near  the  field  boundaries  and  for  email  values  of  kr 
that  the  greatest  deviation  between  the  results  given  by 
the  asymptotic  series  and  the  eigenfunction  series  is 
expected.  In  the  examples,  the  amplitude  and  phase  of 
G-j  were  calculated  by  r.eans  of  the  eigenfunction  series, 
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8? 

and  these  values  wore  taken  to  bo  the  reference  or  true 
ones*  These  calculations  were  performed  on  the 
OSU  -  IBM  -  709^  digital  computer.  In  the  examples 
involving  v/cdges  with  oi  greater  than  180  degrees,  the 
fields  were  also  calculated  by  means  of  the  ordinary 
Pauli  series,  the  generalized  Pauli  ceries,  and  the 
Obcrhe ttinger  series.  In  the  examples  involving  wedges 
with  less  than  180  degrees,  the  fields  were  calculated 
using  the  generalized  Pauli  series  and  the  generalized 
Oberhettingcr  series.  The  results  of  these  computations 
are  compared  v/ith  those  given  by  the  eigenfunction  scries. 
In  all  cases,  the  values  of  have  been  divided  by  2  in 
order  to  bring  them  into  correspondence  v/ith  the  values 
published  by  v.’ait. 

Each  example  is  preceded  by  a  description  of  the 
assumed  geometry  and  brief  summary  of  the  results.  The 
wedge  angle  and  the  angles  of  incidence  and  observation  are 
denoted  on  the  graphs  as  are  the  boundary  conditions.  The 
Neumann  boundary  condition  is  denoted  by  the  acoustically 
equivalent  tern  "hard,"  v/hile  Inc  Dirichlet  condition  is 
described  by  "soft.'1  The  particular  asymptotic  series 
used  in  each  calculation  is  also  denoted  on  the  graphs. 
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The  eigenfunction  result s  are  represented  by  a  heavy  dashed 
lino.  The  asymptotic  scries  results  arc  denoted  by  colic! 
or  dashed  lines  or  by  circles.  The  results  obtained  using 
the  leading  term  as  well  nc  higher  order  terms  of  the  asymp¬ 
totic  series  are  illustrated  on  the  graphs,  The  results 
given  by  the  first  term  of  a  series  arc  denoted  by  (l); 
those  obtained  with  the  first  and  second  by  (2);  by  the 
first,  second,  and  third  by  (2);  cj_c. 

The  results  obtained  for  from  the  "first  tern”  of 
any  of  the  asymptotic  series  include  all  the  plane  wave 
contributions  plus  the  sum  of  two  or  four  Fresnel  integral 
terms  which  describe  the  leading  tern  of  the  particular 
asymptotic  scries.  The  ordinary  Pauli  scries  requires  two 
such  integrals,  while  the  other  scries  require  four.  The 


results  obtained  from  the  "first  and  second  terms"  include 


the  contributions  listed  above  plus  t lie  sum  of  two  or  four 
terms  given  by  the  second  term  in  the  appropriate  scries. 
The  uncorrcctcd  asymptotic  series,  Equation  l*+8,  was  not 
used  in  these  calculations. 


Example  1_ 

In  this  example  a  plane  wave  3s  assumed  to  be  incident 
on  a  hard  w edge  having  an  exterior  angle  220°*  The 
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assumed  angle  of  incidence  is  *  =  0°#  The  field  is 
calculated  along  the  shadow  boundary  at  ^  =  l80°.  The 
first  term  of  the  generalised  Pauli  scries  and  the  first 
two  terms  of  the  ordinary  Pauli  scries  describe  this  field 
accurately  for  values  of  kr  as  snail  as  1.0.  Four  terms 
are  required  in  the  Oberhct tingcr  series  to  produce  good 
agreement  with  the  eigenfunction  series  for  5  ^  kr  <  15» 
BcloV/  kr  a  5  even  four  terns  do  not  produce  good  agreement. 
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Example  J 


In  this  example  the  exterior  wedge  angle  ic  assumed  to’ 
be  190  degree';  the  angle  of  incidence  0  degrees*,  and  the 
field  is  calculated  on  the  shadow  boundary  at  loO  degrees. 

The  Pauli  series  fails  to  yield  an  accurate  value  of  G 

r 

on  the  shadow  boundary  regardless  of  the  number  of  terms 
used.  Two  terms  are  required  in  the  generalized  Pauli  series 
and  the  Oberhettinger  scries  to  produce  good  agreement 
v/ith  the  values  given  by  the  eigenfunction  series  for 
kr  <  5*  When  kr  is  greater  than  5  only  one  terra  is  required 
in  the  generalized  Pauli  series,  but  two  are  required  in 
the  Oberhettinger  series  up  to  kr  =  15» 
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Kx.'ir.nlo  ^ 

The  exterior  wedge  angle  i n  this  cxanpl c.  is  assumed  to 
be  ICO  decree::.  The  angle  of  .incidence  is  10  decrees.  The 
field  is  calculated  at  155  decrees,  midway  betv/een  one  of 
the  facer;  of  the  wedge  and  a  reflection  boundary.  This  is 
the  first  example  involving  a  wedge  of  angle  less  than 
l8j  degrees.  The  field  is  calculated  using  the  eigen func t ion 
eerier;,  generalised  Pauli  series,  and  the  generalised 
Oborho ttingcr  scries.  lioth  hard  and  soft  wedges  are 
examined.  The  generalised  Pauli  and  the  generalised 
Cbcrhct linger  series  are  equally  useful  in  this  case.  The 
generalised  Pauli  is  slightly  belter  since  it  describes 
the  amplitude  of  the  field  mere  accurately  at  small  values 
of  kr  than  the  Cberho ttingcr  series  does.  The  phase  of 

Gd  is  not  shown  in  this  or  the  following  examples. 

* 

It  can  be  shown  that  if  «  =  ,  where  M  is  an  integer, 

the  diffraction  components  of  Gp  vanish  and  the  total  field 
is  described  completely  in  terns  of  plane  wave  components. 
This  corresponds  to  the  fact  that  for  these  particular 
wedge  angles  the  problem  can  be  analysed  in  terms  of  ordinary 


image  theory. 
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Example  £ 

The  assumed  wedge  angle  is  100  degrees  in  this  example; 
the  angle  of  incidence  is  30  decrees;  and  the  field  is 
calculated  alone  the  reflection  boundary  at  ~  10°.  The 
first  tern  of  the  generalised  Pauli  series  accurately  dcccrib 
the  field  for  hr  as  snail  os  1  in  the  case  of  tli c  hard  wedge. 
In  the  case  of  the  soft  vve dee,  two  terns  of  tin's  series  are 
required  to  produce  good  agreemen t  with  the  eigen func tion 
series  at  snail  values  of  hr.  Two  terns  are  required  in 
the  Oberhet linger  series  to  describe  the  field  of  the  hard 
v/edge,  In  the  ease  of  the  soft  wedge,  one  tern  suffices 
except  at  snail  values  of  hr. 
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The  exterior  wedge  angle i  o  nnsu*»iO  d  to  be  70  degrees; 
and  the  angle  of  incidence  10  degrees.  In  this  example,  the 
field  was  calculated  at  ft  ~  60°,  Midway  between  one  face  of 
the  wedge  and  the  reflection  boundary  at  $0°.  In  the 
ease  of  the  hard  v.edgc,  two  tern:;  arc  required  in  the  gen¬ 
eralised  Pauli  series  to  produce  good  agreement  with  the 
eigenfunction  scries  for  hr  <  2.^.  In  the  case  of  the  soft 
wedge,  only  one  tern  is  required.  The  first  tern  of  the 
Oberhet tinge r  series  describes  the  field  of  the  soft  wedge 
quite  well.  This  series  does  not  yield  particularly  good 


results  in  the  ease  of  the  hard  wedge* 
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Example  £ 

In  this  final  example,  the  exterior  wedge  angle  io 
aosumed  to  be  degrees;  the  angle  of  incidence  20  degrees; 
and  the  angle  of  observation  10  degrees*  In  the  case  of  the 
hard  wedge,  the  values  for  Gp  obtained  with  the  first  term 
of  the  generalised  Pauli  and  Cbcrhet tinger  ccricG  are  in 
excellent  agreement  v/ith  those  given  by  the  eigenfunction 
series*  In  this  case,  the  diffracted  component  of  the  field 
is  very  or.nll,  and  a  simple  superposition  of  the  plane  wave 
terns  of  Equation  155  in  sufficient  to  describe  the  field. 

In  the  ease  of  the  soft  v/edge,  the  first  tern  of  the  gen¬ 
eralized  Pauli  series  describes  the  field  accurately  for 
values  of  k.r  greater  than  3*5*  For  smaller  values  of  fcr 
thin  series  fails.  The  generalized  Oberhc ttinger  series 
does  not  describe  the  field  adequately  for  the  soft  wedge 
ease  over  the  range  tested*  This  example  illustrates  that 
both  the  generalized  Pauli  and  the  generalized  Cterhet tingcr 
series  fail  for  snail  values  ofc<,  but  that  the  generalized 
Pauli  series  fails  more  gracefully  than  does  the  generalized 
Oberhc ttinger* 
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Figure  46,  Magnitude  of  Cp  for  a  40  degree,  soft  wedge  (a  comparison  of  the 
generalised  Pauli  and  the  eigenfunction  series).  The  circles  denote  the  values 
given  by  the  first  term  of  the  generalised  Pauli  series.  The  second  and  third 
terns  produce  insignificant  variations  in  these  values. 
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Conclusions 

Tho  important  rcoulto  of  this  work  can  bo  ounr.nrir.od 

mo  follov/o: 

•  i 

•)  The  intocr.il  expression  for  the  total  field  produced 
by  tho  diffraction  of  a  piano  wave  by  a  t\«. -dimensional 
wedge  can  bo  expanded  ocynptoticnlly  by  a  variety  of 
techniques'  but  the  resulting  expansions  can  be  shown 
to  bo  intimately  related*  In  particular,  the  general i nod 
Pauli  ccrlco  and  the  generalised  Oberhettinger  corico 
•re  einply  different  nrrangenente  of  the  cane  series, 
and  the  ordinary  Pauli  ccricc  reprcr.cntc  a  special  form 
of  the  generalised  Pauli  corico* 

b)  Of  the  asymptotic  c erica  examined  in  this  work,  the 
generalised  Pauli  ecrics  io  the  moot  satisfactory  for 
eoaputing  the  field  near  transition  regions*  The  firot 
terra  in  this  series  provides  an  accurate  description 
of  the  amplitude  and  phase  ol  the  field  over  a  wide 
range  of  valuer,  of  «C  and  kr*  It  io  particular.,/ 
ouperior  to  the  Oberhettinger  series  in  cnees  involving 
wedges  having  angles  quite  different  fron  ISO  degrees* 

c)  The  ordinary  Pauli  series  is  ouperior  to  the  Oberhettinger 
corics  for  calculating  the  fields  diffracted  by  wedges 


hnvinc  lornc  exterior  an^len,  while  the  Obcrhcttinccr 
oeriea  in  r.uperior  to  the  ordinary  Pauli  series  for 
ensen  involving  wedges  v/itli  exterior  angles  in  the 
neighborhood  of  and  3cnn  than  180  degrees. 


CHAPTER  V 


DIFFRACTION  COEFFICIENTS 


In  this  chapter  the  leading  terms  of  the  asymptotic  expansions 
given  in  the  summary  of  equations,  pp.  79-85,  are  used  to  determine 
the  diffraction  coefficient  for  the  wedge.  In  terms  of  the  geometrical 
theory  of  diffraction  the  diffracted  electric  field  for  plane  and 
cylindrical  waves  normally  incident  on  the  wedge  is  given  by 

E^r,?)  -  E1  .  [zz.  D  (M'>  -*'i  nh (*,*•)]  ,  (159) 

\lr 


where  E  is  the  electric  field  incident  on  the  edge, 

z  is  a  unit  vector  parallel  tc  the  edge, 

D  ,D,  are  the  scalar  diffraction  coefficients  for  the 
s  n 

Dirichlet  (soft)  and  Neumann  (hard)  boundary 
conditions,  respectively, 

4,0*  and  r  arc  defined  in  Figure  2. 

In  the  discussion  to  follow  E^j  and  correspond  to  the  soft  and 
hard  diffracted  components,  respectively,  of  Gp  in  Equation  147,  pro¬ 
vided  that  E*  and  E^t,  respectively,  equal  one.  At  thin  point  it  is 
convenient  to  introduce  the  definitions 


D'«+*') 


n  JTfry 


sin  7T/n 


4t  0' 


COS  —  -  C08 

n 


n 


+ 


_c  (  7t  +  (0  +  ^')\ 

2n/m  A  5  J 


(160) 


(161) 
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where  the  superscript  +  is  directly  associated  with  the  first  + 
sign  in  the  argument  of  the  cotangent, 


F|kr  a"  (0+  0’)] 


2jjkr  a”  <0  +  V)  c:)Vr  a  >J 


00 


e"lC  d<  (162) 


/  s — 

J  kra 


in  which 


+  + 
a  (0+  (/’)  *  1  +  cos (2rmM  -  («  +  <^')), 


(163) 


where  N  is  the  integer  which  most  nearly  satisfies  the  equation 

+ 

-  (tf +  <*')  +  27mir  =  +7T  .  (164) 


The  superscript  +  is  directly  associated  with  the  +  sign  preceding  the 
77.  The  positive  branches  of  the  square  roots  are  taken  in  Equations 
160  through  162. 

+ 

The  magnitude  and  phase  of  F(kra  )  are  plotted  in  Figure  48, 

+  + 

where  it  is  seen  that  as  kra  increases,  iFl  -V  1.  For  kra  >  10, 

F  fifl. 

From  the  preceding  discussion  and  the  leading  term  in  the  Generalized 
Pauli  asymptotic  series,  Equation  153,  the  set i  >v  diffraction  coefficients 
(*,*•)  “ 

{d+(0-ft')F[kra+((/- 0' )  ]+d‘ (0-/-' ) F[kra’  (0-?’)  ]j 
?^d+(0+0’)K[kra+(0+f)]  +  d*(0+(")F[kra"(f.+0’>]j 


(lf.5) 


134 


in  which  the  upper  sign  between  the  two  bracketed  terms  is  associated 

with  the  soft  (s)  boundary  condition  and  the  lower  sign,  with  the  hard  (h) 

boundary  condition.  At  points  well  removed  from  the  transition  regions 

adjacent  to  incident  and  reflected  field  shadow  boundaries  referred  to 

+ 

as  field  boundaries  earlier,  where  kra  >  10,  the  correction  factors 
F  may  be  replaced  by  unity,  and  Equation  165  simplifies  to 


D 

it 


(*>> 


(166) 


which  is  the  well-known  form  of  the  wedge  diffraction  coefficient  given 
by  Keller."* 

Next,  using  the  leading  term  in  the  Pauli  asymptotic  series, 
Equation  157,  the  scalar  diffraction  coefficient 

D  <M'>  - 
fl 

D’(M')  F[kra(M')l  +  D'(ft>)  F[kraM>+^)];  (167) 

+ 

in  the  Pauli  solution  N  is  zero  so  that 
+ 

»■($+#’)  -  o(4+<h  -  1  +  costy+f).  (168) 

Outside  of  the  transition  regions,  where  kra  >  10,  the  correction 
factors  may  be  replaced  by  unity  and  Equation  167  simplifies  to 
Equation  166. 

Finally  using  the  leading  term  in  the  generalized  Oberhettinger 
series,  Equation  158,  the  scalar  diffraction  coefficients 
0(0,$')  - 

R  -j- 

[d  +(<M')  +  1 

2  J  7Tka*0-  ^  ) 


{l-F[ka+(0-^’)]j 
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,  (169) 


+  + 

where  N  is  determined  from  Equation  164  nnd  a  (y+f)  from  Equation  163. 

+ 

Outside  of  the  transition  regions,  where  kra  >  10,  the  correction 
factors  again  may  be  replaced  by  unity  and  it  is  seen  that  Equation  169 
simplifies  to  Equation  166. 

Three  expressions  have  been  given  for  the  scalar  diffraction  coeffi¬ 
cients,  Equations  165,  167,  and  169.  Unlike  the  earlier  expression  given 
3 

by  Keller  these  are  not  restricted  to  the  region  outside  the  transition 
regions.  It  is  natural  to  ask  which  is  the  most  useful  in  the  transition 
regions,  in  terms  of  its  accuracy  and  simplicity.  It  is  evident  that 
Equation  169  is  the  most  complicated  representation  of  the  diffraction 
coefficient,  whereas  Equation  167  is  the  most  simple.  On  the  other 
hand,  a  study  of  the  examples  in  Chapter  IV  reveals  that  the  most 
accurate  results  (assuming  only  the  leading  term  is  retained  in  each 
asymptotic  series)  are  obtained  with  the  Generalized  Pauli  series. 

A  computer  subroutine  has  been  written  for  the  diffraction  coefficient 
described  by  Equation  165,  and  it  has  been  used  in  a  number  of  wedge 
diffraction  problems^^»123,124  vith  good  accuracy,  provided 
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r  >  0.2  wavelength.  It  is  definitely  more  accurate  than  the  expression 
for  the  diffraction  coefficient  given  in  Equation  167.  However,  a 
computer  subroutine  is  not  essential  to  the  use  of  the  diffraction 
coefficients  based  on  the  Generalized  Pauli  series,  since  the  cor- 
rections  in  the  transition  regions  can  be  easily  made  using  Figure  48. 

The  dyadic  diffraction  coefficient 

d(M’>  -  «v^’>  - 

in  Equation  159  can  be  applied  to  three-dimensional  problems  provided 

that  the  incident  field  at  the  point  of  diffraction  is  normal  to  the 

edge.  Furthermore,  it  is  applicable  to  cases  where  either  the  source 

point  or  the  field  point  is  close  to  the  edge.  This  will  be  described 

126 

in  a  later  report  where  the  Generalized  Pauli  solution  has  been 
extended  to  obtain  a  dyadic  diffraction  coefficient  for  waves  obliquely 
incident  on  a  perfectly-conducting  wedge.  This  dyadic  diffraction  coeffi¬ 
cient  can  be  used  in  the  transition  regions  for  plane,  cylindrical  or 
spherical  wave  illumination  of  the  edge. 


Fig.  48.  The  magnitude  and  phase  of  the  correction  factor. 
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